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INTRODUCTION 


Most people who program microcomputers in BASIC are soon 
disheartened by its shortcomings. Firstly, there are some 
applications which can only be sensibly written in assembly 
language. Secondly, BASIC can be slow. This is not to say that BBC 
BASIC is badly written. On the contrary, it runs faster than most 
other dialects of the language. The slowness is simply inherent in 
the language. Thirdly, programs written in BASIC consume a lot of 
memory. For the BBC Micro, this problem becomes acute in the 
graphics modes which leave precious little space for a program. 
BASIC is an interpreter and not a compiler and many of its 
shortcomings are attributable to that fact. 


A compiler verifies the user program (source code) in a separate 
compilation step, generating a machine code version (usually) of 
that program known as the object code. The object code is stripped 
of all comments. It is the object code which is executed at run time. 


An interpreter, however, acts on the source code at run time. It 
handles the code in a line-by-line fashion, checking syntax and 
parsing each line every time it is executed. This increases the 

execution time of the program. It is also wastes space, since the 
source code has to be present in memory. Thus interpreters use 
more memory and run more slowly than compiled object code. 


The solution to the problem appears to be obvious — buy a 
compiler! There are, however, problems associated with compilers 
also. If the compiler is disk or tape based, then it must be loaded 
into memory for the compilation step (though not for execution). 
This limits the size of the source code that can be written. A ROM 
based compiler gets round this problem. Be wary though! Most of 
the compilers for the BBC micro do not generate machine code 
object programs, but rather an intermediate code which can be 
interpreted into machine code at run time. This means that the 
object code will only run on micros equipped with the same 
compiler. Another problem to be considered with compilers is their 
efficiency. A compiler generates a number of machine code 
instructions for each source language statement. The efficiency of a 
compiler can be thought of as the ratio of the minimum number of 
machine code instructions needed to perform a task, to the number 
generated by the compiler from source code written to perform that 
same task. As a general rule, the more friendly the language, the 
less efficient the compiler. 


There is a compiler that all BBC micro owners possess. It's called 
the assembler. The distinction between assemblers and compilers is 
that the former have a one-to-one relationship with the machine 
code that they generate; that is to say each assembly language 
command generates a single machine code instruction. The 
drawback is that assembly language programming can be difficult, 
time-consuming and error prone. Moreover, assemblers do not have 
in-built commands to handle sines, cosines, square roots, random 
numbers etc. 


On the face of it, the user appears to be caught in a cleft stick. The 
user can either choose BASIC for its ease of use, tolerating speed 
and size problems, or assembly language, foregoing floating-point 
arithmetic, trigonometry and so forth. 


Fortunately, the BBC BASIC interpreter consists of a large number 
of small, machine code subroutines, many of which could usefully 
be invoked from an assembly language program. In this book, 69 
such subroutines are described covering 32-bit integer arithmetic, 
floating-point arithmetic, trigonometry and so forth. 


The approach has several advantages: 
а) The subroutines are ROM based and occupy no valuable RAM. 
b) The subroutines are tried and tested. 


c) Functions such as sine, cosine, square root and random 
numbers can be used in assembly language programs without 
the need to write this code from scratch. 


d) The subroutines frequently incorporate useful error reporting. 


e) Object code will run faster and occupy less memory than the 
equivalent BASIC code. 


It is only fair to point out the disadvantages of this approach: 


a) This book covers BASIC 1 and BASIC 2 and the technique 
described will work only on these two ROMs. It will not work 
with HI BASIC (6502 second processor) or with US BASIC (US 
BBC micros). Neither will it work with any future releases of 
BASIC that may be issued. Consequently, it would be most 
unwise to use the technique directly in any program which is 
destined for sale to the public, or where the user plans to 
upgrade the BASIC ROM within the life time of the program. 
Even in these cases, however, this book has distinct value. 


Many of the subroutines іп the BASIC ROMs are written with 
an elegance and tightness which is unlikely to be surpassed by 
the user. A study of these subroutines will undoubtedly profit 
the user who is forced to write similar code. There will be 
many instances, however, when a user may wish to use the 
technique directly. It would be a pity to deprive all users of an 
effective and time-saving technique, simply because the 
technique is not suitable for every occasion. 


b) Because assembly language is used, the source code is often 
more long-winded to write than the equivalent BASIC. 


с) The technique is not a cure-all for all problems. Although the 
user is offered an easy path to sines and cosines, not much 
more than 1096 of the execution time can be saved for these 
complicated functions. Chapter 9 provides some useful tips in 
this respect. 


d) The user must have some acquaintance with assembly language 
programming and such knowledge is assumed in this book. 


This then is not a primer in assembler. The bulk of this book 
contains descriptions of BASIC subroutines, their entry points, 
timings and set-up conditions. This is supported by fully 
documented and tested examples of code, making it possible for 
relatively inexperienced assembly language programmers to use the 
techniques to develop quite sophisticated applications. 


The book also contains in places descriptions of the theory 
necessary for a full understanding of the technique advocated. АП 
programmers inevitably make mistakes and this knowledge is 
indispensable when debugging those mistakes. 


Much of this book is to do with numbers and their representation 
within the BBC micro. A feel for binary and hexadecimal 
numbering systems is necessary to understand completely the text 
presented here. In particular, floating point numbers require that 
the concept of the binary point be grasped. Therefore, the next 
chapter is devoted to numbering systems. Experienced 
programmers will doubtless skip over this chapter. It is specifically 
intended for the less experienced programmer and for this reason 
emphasises the ‘why’ as well as the ‘how’ of numbering systems. 


1 NUMBERING SYSTEMS 


Since it is in every day use, the decimal numbering system is 
widely understood. It is based on (has a radix of) ten. The radix of a 
numbering system embodies several important concepts: 


a) itis the number of possible values that a single digit can 
contain. (decimal values are 0 to 9, ten values in all). 


b) itis one more than the maximum value of a single digit. 


c) each digit in the number will represent some power of the 
radix. It is the position of the digit within the number which 
determines the power of the radix by which the digit is 
multiplied. For example, the decimal number 1932.74 - 


10*10*10 10x10 10 1 1/10 1/100 
1 9 3 2 7 4 


It will be shown that the decimal system is not ideally suited to 
computers. Fortunately, the alternative numbering systems which 
suit computers best use the same principles as the decimal system 
and so they are not difficult to master. 


1.1 The Binary System 


A computer is an electronic machine. It knows nothing of numbers. 
It consists of many components, each of which is designed to 
respond to pulses of electricity, providing the pulses arrive in a 
particular pattern and at a particular time. At the heart of the 
computer lies a central processor unit (cpu) which not only 
performs the arithmetic functions requested by a program, but also 
exercises control over other component parts of the computer. 


At any given instant in time, an electrical pulse in a circuit may 
either be present or absent. In other words, a circuit can be in one of 
two states. This is a binary (meaning based on two) system. It is 
human beings that allocate numbers to these two states. The binary 
system has only two numbers, 0 and 1. 0 is taken to mean that the 
pulse is absent, whilst 1 denotes its presence. 


Ав ап example of these pulses, each cpu has а repertoire of 
commands that it can perform, known as its instruction set. At 
specific intervals of time, the cpu will interpret a pattern of pulses 
as a command to perform a particular function. This function could 
be to add or to store etc. The pulses arrive simultaneously down 
eight wires which collectively are called the data bus. For example: 


data bus pulse binary 
wire 7 present 1 this 
wire 6 absent 0 pattern 
wire 5 absent 0 оғ 
міге 4 absent 0 pulses 
wire 5 present 1 is 
wire 2 absent 0 the 
wire 1 absent 0  DEY 
wire 0 absent 0 instruction 


A computer program causes pulses just like this to move from one 
part of the computer to another (albeit with a good deal of help 
from both language software and the operating system). In the 
example above, the individual pulses are known as bits (binary 
digits) and the collection of 8 bits that were sent down the data bus 
is called a byte. 


Thus a computer, its design constrained by the laws of physics, 
operates in a binary fashion. It follows that it will be easier to work 
with computers if the programmer achieves some expertise with the 
binary system also. 


In the binary system, each digit represents some power of two, such 
that the binary number 01101000, for example, is equivalent to: 


(0x128) + 


(1x64) + 
(1x32) + 
(0x16) + 
(1x8) + 
(0x4) + 
(0x2) + 
(0x1) 


= 104 in decimal. 
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1.1.1 Binary Addition 


Binary addition follows rules analogous to decimal addition. Each 
time the sum of a column is two or more (rather than ten or more) 
there is a carry to the next column. Since there are only two digits 
in the binary system, it is possible to define very simple rules for 
binary addition: 


00 =0 
01 -1 
1+1 =O (carry 1) 


1+1+1 = 1 (carry 1) 

These simple rules can be applied to much larger binary numbers. 
For example: 

128 64 32 168 4 21 decimal 


01101000 104 
01011101 93 + 
11000101 197 


1.1.2 Binary Subtraction 
Simple rules can also be defined for binary subtraction: 
0-0 - 0 


1-0 - 1 
0-1 = 1 (borrow 1) 
1-1 = 0 


Likewise these simple rules can be applied to the subtraction of 
larger binary numbers. For example: 


128 64 32 168 4 21 decimal 


01101000 104 
01011101 95 - 
00001011 11 
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1.1.3 Negative Binary Numbers 


So far only positive, binary integers have been considered. Negative 
binary integers are held in twos complement form, for reasons 
which will soon become apparent. To convert a binary number to 
twos complement form, it is only necessary to change all zeroes to 
ones and all ones to zeroes, adding 1 to the result. In the example 
above, the binary for decimal 93 was seen to be 01011101. 


To obtain the binary for –93: 
01011101 = +93 


10100010 = change 0 to 1 and vice-versa 
00000001 add 1 


= 
о 
= 
e 
e 
e 
= 
= 
Ш 


-93 


1.1.4 Binary Subtraction By Addition 


The twos complement form for negative numbers makes 
unnecessary any subtraction circuitry in the cpu. Instead, 
subtraction can be achieved via the cpu's adder. In the example 
used for subtraction, 104 - 93 is the same as 104 + (-93). Therefore 
by converting 93 to its twos complement form, addition can be used 
to achieve subtraction providing any final carry is ignored: 


01101000 104 
10100011 -93 + 


ignore carry 1 00001011 11 


1.1.5 Binary Fractions 


Binary fractions follow concepts analogous to decimal fractions. In 
decimal fractions each digit position represents some power of 
1/10. Thus the decimal fraction 0.875 = 


8 7 5 


— + — + —— 


10 100 1000 


With binary fractions, each digit position represents some power of 
1/2 and the point is known as the binary point (rather than decimal 
point) to emphasise this fact. Thus the binary equivalent of 0.875 is 
0.111 = 


12 


1.2 Тһе Hexadecimal System 


Binary numbers have ап obvious disadvantage. It takes а large 
number of zeroes and ones to represent even comparatively small 
decimal numbers. But computers work in a binary fashion and 
conversion between the two radices is a laborious process. What is 
required is some shorthand version of binary. 


For a numbering system to relate directly to binary, it must be based 
on some power of two. The choice of numbering system is mainly 
determined by the number of bits in a byte. In practice the radix 
chosen fulfils the following conditions: 


a) itis a power of 2 

b) itisas close to 10 as possible 

c) itsubdivides a byte into equal proportions (usually halves but 
not always). 


A BBC micro has 8 bits in a byte. Therefore, a byte subdivides into 
two equal size quartets. This fixes the numbering system as 
hexadecimal with a radix of 16 (one more than the maximum 
decimal number that can be stored in four bits). Had there been 6 
bits in a byte, the byte would have been regarded as two triplets and 
the numbering system would have been octal (based on eight). 


There are sixteen digits in the hexadecimal system and a symbol is 
required to represent each digit. Clearly for the digits 0 to 9, the 
same symbols (0 to 9) can be used. The single hexadecimal digits 
that represent decimal numbers 10 to 15 are a bit more of a 
problem. In fact the symbols A to F are allocated to these numbers. 
When used in this way, these symbols should not be confused with 
the ASCII letters ‘A’ to ‘F’. The meaning is always understood 
because in both BASIC and assembly language, hexadecimal 
numbers are preceded by an ampersand (&). 


Thus in a single byte, the range of hexadecimal numbers that can be 
stored is from &00 to &FF (0 to 255), with each quartet holding from 
%0 to «ХЕ (0 to 15) as shown below: 


binary hex dec binary hex dec 
0000 %0 0 1000 %8 8 
0001 &1 1 1001 %9 9 
0010 %2 2 1010 &A 10 
0011 &5 5 1011 &B 11 
0100 %4 4 1100 &C 12 
0101 %5 5 1101 @ 13 
0110 %6 6 1110 &E 14 
0111 %7 r4 1111 &F 15 
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1.2.1 Binary/Hexadecimal Conversion 


The most important property of a hexadecimal number is that it 
can be derived instantly from a binary number. To do this, the 
binary number is subdivided into quartets starting at the least 
significant end. If the most significant bit(s) are not a quartet, they 
should be padded with leading zeroes to make an exact quartet. 
Hexadecimal can then be substituted for each quartet individually. 
For example, consider the binary number 11101011010010011. 
Splitting into quartets gives: 

0001 1101 0110 1001 0011 
hex = 1 D 6 9 3 


= &1D693 


The converse operation from hexadecimal to binary is equally 
simple. It is impossible to do anything во simple with decimal 
numbers. 


1.2.2 The Roundness of Hexadecimal 


There is another bonus obtained by using hexadecimal. Because the 
design of the computer is based on the number 2, hexadecimal 
frequently results in an easily remembered, round number where 
the decimal equivalent does not. One example of this is found in 
memory addressing: 


total BBC model В memory 10000 65556 


1K of memory 400 1024 
one page of memory 100 256 
PAGE (по disks/Econet) Е00 3584 
start of BASIC 8000 32768 
start of MOS с000 49152 


Another example сап be found іп the ASCII character set: 


From &00 to &1F = control characters CVDU O to VDU 31) 
From 850 to 839 = 0 to 9 
From &41 to 85A = A to Z 
From 861 to &8A = а to 2 


It will be seen that an ASCII number is the number plus &30. The 
upper-case letters are the number of the letter in the alphabet plus 
&40, whilst lower-case letters are the number of the letter in the 
alphabet plus &60. The control characters, representing VDU 0 to 
VDU 31, can be entered from the keyboard by using CTRL and 
another key. The action of the CTRL key is to subtract &40 from 
that other key. Thus VDU 1 = CTRL A, VDU 2 = CTRL B etc. 
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1.2.3 Hexadecimal Addition 


In hexadecimal addition there is a carry to the next more 
significant digit whenever the sum exceeds decimal 16. Thus the 
highest number that can exist in any digit is &F (15). Consider the 
addition of &A and &B, for example. Since &A is equivalent to 
decimal 10 and &B is equivalent to decimal 11, the decimal sum is 
21. As this is greater than decimal 16, 16 must be subtracted and a 
carry must be generated to the next more significant digit. Thus the 
hexadecimal sum is &15. The right hand digit (5) is obtained from 
21 – 16, whilst the left hand digit is the carry. 


1.2.4 Negative Hexadecimal Numbers 


The hexadecimal complement form is derived by subtracting the 
number from a string of &F's and then adding 1. The number of 
&F's in the string depends on the precision of arithmetic in 
question. Thus in 16 bit arithmetic, &FFFF would be used, whilst 
in 32 bit arithmetic &FFFFFFFF would be used. As an example, the 
16 bit hexadecimal complement of &ABC is derived as follows: 


1.2.5 Hexadecimal Fractions 


In hexadecimal fractions, each digit position represents some 
power of 1/16. The point is known as the hexadecimal point (rather 
than decimal point) to emphasise this fact. Thus the decimal 
fraction 0.875, which is 14/16, has a hexadecimal equivalent of 0.E. 


15 


16 


2 INTEGERS 


In BASIC АП integer fields occupy four bytes (32 bits). This allows 
for positive integers from %00000000 to &7FFFFFFF (0 to 
2,147,483,647). Negative numbers range from &FFFFFFFF to 
&80000000 (-1 to -2,147,483,648). 


BASIC recognises ап integer variable by the % at the end ofits 
name. Some of these variables, known as the resident integer 
variables, @% and A% to Z%, occupy fixed locations in RAM. 
Other integer variables are located as referenced in an area of RAM 
following the program text. The assembler programmer is free to 
use the resident integer variables. As with BASIC, these memory 
areas can be used to pass parameters from one program to another. 
However, 0% and P% have special significance as location 
counters in assembly language and must not be used (0% can 
safely be used by BASIC 1 users). It is also inadvisable to use @% 
which controls BASIC print formatting. The assembler program 
accesses the resident integer variables by address rather than by 
name. А list of these addresses may be found in Chapter 7. 


2.1 Integer Work Areas 


Integer numbers are always stored with the least significant byte 
first and the most significant byte last. Thus the number &12345678 
would be held as: 


byte 0 = 878 
byte 1 = 856 
byte 2 = 834 
byte 3 = 812 


The BASIC interpreter performs all of its integer arithmetic in four 
bytes of working storage in Page Zero. The four bytes are &2A, &2B, 
&2C and &2D. From henceforth, these four bytes will be referred to 
as the Integer Working Area or IWA. These four bytes are also used 
by the interpreter for other purposes. When used as the IWA, the 
normal rules for integer variables are obeyed. &2A contains the 
least significant byte, whilst &2D contains the most significant byte. 


BASIC has its own stack located immediately below HIMEM. Like 
the processor stack, it runs downwards in memory. The stack is 
maintained by a set of subroutines within the BASIC ROM. A stack 
pointer is held in &4 and %5 (lo,hi). Generally, whenever BASIC 
has two integer fields to process, one will be located in the IWA 
while the other is in the stack, pointed to by &4 and &5. When 
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using BASIC routines we shall usually load one integer into the 
IWA and point &4,&5 at the other. 


BASIC also has to keep track of whether it is processing integers, 
floating points or strings. It achieves this in two ways. Firstly 
subroutines return a value in the A register as follows: 


А = 0 processing a string 
A = %40 processing an integer 
А = &FF processing floating point 


It also stores the type of variable in &27 as follows: 


&27 = 0 byte 

&27 = 4 4 byte integer 

&27 = 5 5 byte floating point 
827 = &81 string 

&27 = &A4 function 


&27 = &F2 procedure 


From time to time, when using the 32-bit integer subroutines,it 
will be necessary to set either registers or memory areas to 
conform with the above. 


2.2 Defining Integer Constants 


For BASIC 2 owners, defining a 32-bit integer constant to the 
assembler is a matter of the utmost simplicity. The assembler 
directive EQUD is provided for this purpose. For example: 


10 DIM mc% 100 

20 FOR pass% = 0 TO 2 STEP 2 
30 PX = mc% 

40 E 

50 ОРТ pass% 

60 .constant EQUD 5000 

70 1 

80 NEXT pass% 

90 STOP 


EQUD is an assembler directive. This means that it instructs the 
assembler to perform a task at assembly time. This is quite different 
from an instruction mnemonic. Mnemonics are translated into 
machine code for execution at run time. 


The EQUD directive instructs the assembler to reserve four bytes of 
memory at the current value of the location counter (P%). It 
automatically reverses the storage of data such that the least 
significant byte is stored first. Thus, since 5000 is &1388, the 
directive stores &88 at the current value of the location counter, 
&13 at the next address, and &00 at the next two addresses. It also 
steps the location counter by 4. 
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BASIC 1 users have to improvise to achieve the same effect. Since 
in BASIC 1 there is only one directive available, namely OPT, this 
must be pressed into service to provide an equivalent mechanism. 
Consider the following: 


10 DIM mc% 100 

20 FOR pass% = 0 TO 2 STEP 2 

30 PZ = тс# 

40 E 

50 ОРТ pass% 

60 .constant \ Label 
70 ОРТ FNEQUD(5000) \ pseudo-directive call 
80 1 

90 NEXT pass 

100 STOP 

110 DEF FNEQUDCAZ) 

120 1P% = AX 

130 PZ = Р + 4 

140 = passX 


At line 70, OPT is made to call a BASIC function. The value of OPT 
must not be changed, so the function must return the current value 
of ОРТ and indeed line 140 does this. The ОРТ directive ensures 
that the call to the function EQUD is performed at assembly time 
(not execution time). This is not a hybrid program. Once 
assembled, the BASIC function is no longer required. Within the 
function itself, statements may be included as required. In the 
EQUD function, the required constant is stored at the current value 
of the location counter. Because an integer variable, A96, is used to 
hold the argument (in this example 5000) passed by the assembler, 
the constant is automatically aligned with the least significant byte 
first. The counter is then stepped by 4. The constant thus set up 
may be referenced by the label on line 60. 


The function, FNEQUD, is an example of a pseudo-directive, a 
technique extensively used in this book. BASIC 2 also has the 
following assembler directives: 


EQUB for 1 byte 
EQUW for 2 bytes (least significant byte first) 
EQUS for strings 


BASIC 1 users can use the equivalent pseudo-directive functions 
below. Note that if the constant supplied to the function is too big, 
it will be truncated at the most significant end. 


Both BASIC 1 and BASIC 2 users will benefit from the RESB 

pseudo directive which reserves a specified number of bytes (1st 

argument) and fills them with a specified character (2nd argument). 
10 DIM mc% 100 


20 FOR pass% = 0 TO 2 STEP 2 
30 P% = mc% 
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40 E 

50 ОРТ pass% 

60 .constant1 

70 ОРТ FNEQUBC50) 

80 .constant2 

90 ОРТ FNEQUW(500) 

100 .string1 

110 ОРТ FNEQUSC"Example of a string") 
120 .reserve 

130 OPT FNRESB(20,0) \ reserve 20 bytes of zero 
140 1 

150 NEXT pass% 

160 STOP 

170 DEF FNEQUBCAZ) 

180 ?PZ = A% 

190 РА = РА +1 
200 = pass% 
210 DEF FNEQUWC(AZ) 


220 ЭРХ = А% МОр 256 
230 ?(P%+1) = AZ DIV 256 
250 P% = PZ + 2 


250 = pass% 

260 DEF FNEQUS(A$) 

270 $PZ = AS 

280 P% = PX + LEN(A$) 
290 = pass% 

300 DEF FNRESB(A%,B%) 
310 FOR I% = 1 TO AX 
320 2P% = B% 

330 PZ = PX +1 

340 NEXT 

350 = pass% 
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2.3 Integer Routines Summary 


The following table summarises 32-bit integer routines available 
within the BASIC ROM. Conversion routines, e.g. integer to ASCII, 
are the subject of a later chapter. Unlike the floating point routines 
covered in the next chapter, the integer routines do not amount to a 
great deal of code. Moreover, set up procedures are sometimes 
more laborious than for floating point. In some applications, it may 
well be preferable to include the integer code within the user 
program, at the same time changing it slightly to make it more 
specific to the application. In this event, the user is advised to 
study the code available within the BASIC ROM, rather than 
re-invent it. 


Name BASIC 1 BASIC 2 Function 
address address 

icomp ФАрВ5 &AD93 IWA = -ІМА 

idiv &9DE7 &9EOA IWA - IWA DIV integer variable 

iin &B365 &B336 Copy integer variable to IWA 

iminus &9C9D &9CC2 IWA = integer variable - IWA 

imod &9DDE &9E01 IWA = IWA MOD integer variable 

imult &9D4A &9D6D ША = ША * integer variable 

ineg1 &ACEA &ACC4 IWA = -1 

jout &B4F2 &B4C6 Copy IWA to integer variable 

iplus &9C36 &9C5B IWA = IWA + integer variable 

ipos &AD94 &AD71 Make IWA positive 

ismall &AF19 &AEEA IWA = 256*Y + А 

itest &9A85 &9AAD test integer variable = < > IWA 

izero  &AEF9 &AECA IWA = zero 

izpin — &AF85 &AF56 Copy integer variable in zero page 
to IWA 

izpout &BE5C &BE44 Copy IWA to integer variable in 
zero page 
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2.4 Integer Routines Description 


The following pages show how to use each integer routine. Each of 
the routines handles signs correctly. For example, (-3)*(-9) gives 
the answer +27. 


This assumes that the set up procedures specified for each routine 
are followed carefully. To ensure that there is no ambiguity, each 
routine is illustrated by a program. These programs are for 
demonstration purposes only and are not meaningful applications. 
They all use BASIC to print results partly because the material 
necessary to avoid this comes later in the book, and partly to 
present the technique in a way that is most simple to understand. 
АП of the programs are written in BASIC 2, but conversion to 
BASIC 1 involves only: 


a) replacing all routine addresses by their BASIC 1 equivalents. 
b) using pseudo-directives instead of directives. 


The use of most of the BASIC ROM’s integer routines is natural, 
since the routines are structured in a way ideally suited to this 
purpose. However, DIV and MOD are exceptions and the method of 
using them has had to be contrived. 


An important point arises from the interpreter's habit of re-using 
zero page locations for different purposes. This can cause 
difficulties in hybrid BASIC/assembly language programs. The 
following rules should be observed: 


a) do not attempt to call these routines from BASIC language. 


b) in hybrid programs, make sure that the results of any 
calculations are safely stored in memory areas within the 
domain of the program. Do not leave data in BASIC's zero page 
areas, if this data will be required later. 


The following short program illustrates the sort of difficulties 
which can arise: 


10 1%2А - 1000 
20 PRINT !&2A 
30 END 
>RUN 

262186 


These routine addresses, therefore, are provided specifically to 
assist assembly language programming and this is the 
recommended way of using them. 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: icomp 

: IWA = -IWA 
BASIC 1 address : 
: &AD93 


&ADB5 


: IWA contains a 32-bit integer 
: IWA complemented 

: A destroyed 

: X unchanged 

: Y destroyed 

: P destroyed 


: 31 microseconds 
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icomp demonstration (BASIC 2) 


0 icomp = &AD93 

10 DIM mc% 200 

20 FOR pass% = 0 TO 2 STEP 2 

30 PZ = тс# 

40 СОРТ pass% 

50 .constant EQUD -1234 N set up value 


60 .result EQUD 0 \ result 

70 .start 

80 LDX #0 ÓN zeroise Loop counter 
100 .Loop1 


110  LDA constant,X \ get next byte of constant 
120 STA &2A,X \ save in IWA 
130 ІМ \ bump Loop counter 
140  CPX #4 \ end of Loop ? 
150 ВСС Loop1 \ no - back 
N 
\ 


160 JSR icomp complement IWA 


170 LDX #0 zeroise Loop counter 
180 .Loop2 

190 LDA &2A,X \ get next byte of IWA 
200 STA result,X \ save in result 

210  INX \ bump Loop counter 
220  CPX #4 \ end of Loop ? 

250 ВСС Loop2 \ no — back 

240 RTS 

250 1 


260 NEXT pass% 
270 CALL start 
280 PRINT !result 
290 END 


>RUN 
1234 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


error reports 


typical timing 


: idiv 

: IWA = IWA DIV integer variable 
: «ОПЕ? 

: &9E0A 

: IWA contains dividend 

: А% contains divisor 


: &19,&1A (lo,hi) point to a string, 


* A96 -RETURN 


: &1В=0 

: &04,&05 (lo,hi) = НІМЕМ 

: À = #&40 

: DIV and MOD are not easy to extricate from the 


BASIC ROM and the set up procedures are 
necessarily somewhat contrived 


: IWA = quotient 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: division by zero 


: 794 microseconds 
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idiv demonstration (BASIC 2) 


10 idiv = &9EOA : 


DIM mc% 200 


20 FOR pass% = 0 TO 2 STEP 2 
COPT pass% 
-dividend EQUD 26 

EQUD 3 

EQUD &0D254120 
EQUD 0 


30 P% = mc£ : 


40 

50 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 


-divisor 
. fudge 
.result 
-start 


. Loop1 


. Loop2 


LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDX 
STX 
LDA 
STA 
LDA 
STA 
INX 
CPX 
BCC 
LDA 
JSR 
LDX 
LDA 
STA 
INX 
CPX 
BCC 


RTS: 


340 NEXT pass : 


>RUN 


CALL start : 


&6 

&4 

&7 

&5 

#fudge MOD 256 
&19 

#fudge DIV 256 
&1A 

#0 

&1B 

dividend,X 
&2A,X 
divisor,X 
&404,X 


#4 

(оор1 
#&40 
idiv 

#0 

&2A,X 
result,X 


#4 
Loop2 
1 


\ 


ие се не „== = жо „ы mL "е PT ME ll 


" AZ" + RETURN 


get HIMEM Lo 

set &4 

get HIMEM hi 

set &5 

point &19 

at fudge Lo 

point &1A 

at fudge hi 

zeroise loop counter 
clear &1B 

next byte of dividend 
save in IWA 

next byte of divisor 
save in AZ 

bump Loop counter 
end of Loop ? 

no - back 

set integer 

call idiv 

zeroise loop counter 
next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 


PRINT !result : END 
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subroutine name : iin 
function : IWA = integer variable 
BASIC 1 address : &B365 
BASIC 2 address : &B336 
entry conditions : &2A,&2B (lo,hi) points to integer variable 
exit status : IWA set up 
: А destroyed 
: X unchanged 
: Y destroyed 
: P destroyed 
typical timing : 34 microseconds 


special note : Note that iin is shown for completeness only. 
The code represented by iin is trivial. The 
following code can be substituted: 


80 \ CODE TO COPY A 32-BIT INTEGER, 


90 N intvar, INTO THE IWA. 
100 .intvar EQUD 12345678 


110 LDX #0 \ zeroise loop count 
120 .Loop 

130 LDA intvar,X \ next byte of intvar 
140 STA &2A,X \ save in IWA 

150 INX \ bump Loop count 

160 CPX #4 \ end of loop ? 

170 BCC Loop \ no – back 
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їп demonstration (BASIC 2) 


0 iin = &B336 
10 DIM mc% 200 
20 
30 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 1 

220 NEXT pass% 
230 CALL start 
240 PRINT !result 
250 END 


>RUN 


PX = mc% 
COPT pass% 
-constant 
.result 
-start 
LDA 
STA 
LDA 
STA 
JSR 
LDX 
. Loop 
LDA 
STA 
INX 
CPX 
BCC 
RTS 


EQUD 1234 
EQUD 0 


#constant MOD 256 
&2A 
#constant DIV 256 
&2B 
jin 
#0 


&2A,X 
result,X 


#4 
[оор 


1255 


ҒОК ра55% = 0 ТО 2 5ТЕР 2 


ме е” 


и” Т” Z Z Z Z 


CCC 


set up value 
result 


LSB of address of constant 
set up &2A 

MSB of address of constant 
set up &2B 

call iin 

zeroise loop counter 


get next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: iminus 

: IWA = integer variable - IWA 
BASIC 1 address : 
: &9CC2 


&9C9D 


: IWA contains a 32 bit integer 

: &04,&05 point to integer variable 
: X = #4 

: IWA set up 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 52 microseconds 
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iminus demonstration (BASIC 2) 


0 iminus = &9CC2 

10 DIM mc% 200 

20 FOR pass% = 0 TO 2 STEP 2 
30 PZ = mc% 

40 СОРТ pass% 


50 .constant EQUD 26 
60 .subtrahend EQUD 3 
70 .result EQUD 0 
80 .start 

90 LDX #0 

100 .Loop1 

110  LDA subtrahend,X 
120 STA &2A,X 

130 INX 

140  CPX #4 

150 ВСС Loop1 

160 LDA #constant MOD 256 
170 STA &4 

180 LDA #constant DIV 256 
190 STA &5 

200  LDX #4 

210  JSR iminus 

220 LDX #0 

230 .Loop2 

240 LDA &2A,X 

250 STA result,X 

260 INX 

270  CPX 44 

280 ВСС Loop2 

290 RTS 

300 1 


310 NEXT pass% 
320 CALL start 
330 PRINT result 
340 END 


>RUN 
23 


A A 


шен тын М те t O uar т е Pd 


CCC 


set up value 
subtrahend 
result 


zeroise loop counter 


get next byte of subtrahend 
save in IWA 

bump Loop counter 

end of Loop ? 

no - back 

LSB of address of constant 
set up 84 

MSB of address of constant 
set up &5 

set up X 

call iminus 

zeroise loop counter 


get next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


error reports 


typical timing 


: imod 

: IWA = IWA MOD integer variable 
: &9DDE 

: &9Е01 

: IWA contains dividend 

: A96 contains divisor 


: &19,&1A (lo,hi) point to a string, 


* A96' + RETURN 


: &1В=0 

: &04,&05 (lo,hi) - HIMEM 

: À = #&40 

: DIV and MOD are not easy to extricate from the 


BASIC ROM and the set up procedures are 
necessarily somewhat contrived 


: IWA = remainder 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: division by zero 


: 782 microseconds 
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imod demonstration (BASIC 2) 


10 imod = &9E01 
20 FOR pass% = 0 TO 2 STEP 2 
30 PX = тсй : 


40 

50 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 


340 NEXT pass : 


>RUN 


.dividend EQUD 26 
.divisor EQUD 3 
. fudge 
.result  EQUD 0 
.start LDA &6 
STA &4 
LDA 87 
STA &5 


LDA Zfudge MOD 256 


STA &19 


LDA Zfudge DIV 256 


STA &1A 
LDX #0 
5ТХ &1B 


. Loop1 LDA dividend,X 


STA &2A,X 


LDA divisor,X 


STA &404,X 
INX 
CPX #4 
BCC Loop1 
LDA Z&40 
JSR imod 
LDX #0 

. Loop2 LDA &2A,X 


STA result,X 


INX 

CPX #4 
BCC Loop2 
RTS:1 


: DIM mc% 200 


COPT pass% 


EQUD &0D254120 


CALL start : 


N 


a tn, mc t aM ыт Sa uit е m қ VI LP өн м” е” uM auo uir 


PRINT 
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" AX" + RETURN 


get HIMEM Lo 

set &4 

get HIMEM hi 

set &5 

point &19 

at fudge Lo 

point &1A 

at fudge hi 
zeroise loop count 
clear &1B 

next of dividend 
save in IWA 

next of divisor 
save in AZ 

bump Loop counter 
end of Loop ? 

no - back 

set integer 

call imod 

zeroise loop count 
next byte of IWA 
save in result 
bump Loop count 
end of Loop ? 

no - back 


lresult : END 


subroutine name 


function 


BASIC 2 address 


entry conditions 


Comments 


exit status 


typical timing 


: imult 

: IWA = integer variable * IWA 
BASIC 1 address : 
: &9D6D 


&9D4A 


: IWA contains a 32 bit integer 
: &04,&05 point to integer variable 
: &27 = #4 


: if the IWA contains an integer larger than 


&FFFF, it is truncated to 16 significant bits. 


: IWA set up 

: A destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 


: 164 microseconds 
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imult demonstration (BASIC 2) 


10 imult = &9D6D : DIM mc% 200 
20 FOR pass% = 0 TO 2 STEP 2 
30 PZ = тс# 

40 СОРТ pass% 

50 .multiplicand EQUD 26 


60 .multiplier EQUD 3 
70 .result EQUD 0 
80 .start 

90 LDX #0 

100 .Loop1 

110 LDA multiplier,X 
120 5ТА %2А,Х 

130 INX 

140 CPX #4 

150 ВСС [оор1 

160 LDA #multiplicand MOD 256 
170 STA &4 

180 LDA Zmultiplicand DIV 256 
190 STA &5 

200  LDX #4 

210  STX &27 

220  JSR imult 

230 LDX #0 

240 .Loop2 

250  LDA &2A,X 

260 STA result,X 

270 INX 

280  CPX #4 

290 ВСС Loop2 

500 RTS 

310 1 


320 NEXT pass% 
330 CALL start 
340 PRINT !result 
350 END 


>RUN 
78 


\ set up value 
multiplier 
result 


—-— 


ет 


zeroise loop counter 


get next byte of multiplier 
save in IWA 

bump Loop counter 
end of Loop ? 

по - back 

LSB of address 

set up 84 

MSB of address 

set up &5 

set up 

&27 

call imult 

zeroise loop counter 


wa a aeo Р Р э э „ы е ше" 


get next byte of IWA 
save in result 

bump Loop counter 
end of [оор ? 

no - back 


CCC 
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subroutine name : ineg1 
function : IWA = -1 
BASIC 1 address : &ACEA 
BASIC 2 address : &ACC4 
entry conditions : none 
exit status : IWA = -1 
: А destroyed 
: X unchanged 
: Ү unchanged 
: P destroyed 


typical timing : 20 microseconds 
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ineg1 demonstration (BASIC 2) 


0 
10 
20 
30 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 


>RUN 


ineg1 = &ACC4 
DIM mc% 200 
FOR pass% = O TO 2 STEP 2 
PZ = тс# 
ГОРТ pass% 
.result EQUD 0 
.start 
JSR ineg1 
LDX #0 
. Loop 
LDA &2A,X 
STA result,X 
INX 
CPX #4 
BCC Loop 
RTS 
1 
NEXT pass% 
CALL start 
PRINT !result 
END 


Ludi" т adi сай 


result 


call ineg1 
zeroise loop counter 


get next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


special note 


80 \ CODE TO COPY 


: iout 

: integer variable = IWA 
: &B4F2 

: &B4C6 

: &37,&38 (lo,hi) points to integer variable 
: &39 must be non-zero 
: IWA set up 

: A destroyed 

: X unchanged 

: Y destroyed 

: P destroyed 

: 37 microseconds 


: Note that iout is shown for completeness only. 


The code represented by iout is trivial. The 
following code can be substituted: 


THE IWA INTO A 


90 N 32-BIT INTEGER VARIABLE, intvar. 


100 .intvar EQUD 0 
110 LDX #0 

120 .Loop 

130  LDA &2A,X 
140 STA intvar,X 
150 INX 

160  CPX #4 

170 ВСС Loop 


\ zeroise loop count 


\ next byte of IWA 
\ save in intvar 

X bump Loop count 
\ end of Loop ? 

\ no - back 
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iout demonstration (BASIC 2) 


0 iout = &B4C6 
10 DIM mc% 200 


20 FOR pass% = 0 TO 2 STEP 2 


30 P% = mc% 
40 СОРТ pass% 
50 .constant 
60 .result 

70 .start 

80 LDX #0 
90 .Loop 

100 LDA constant,X 
110 STA %2А,Х 

120 INX 

130 CPX #4 

140 ВСС Loop 


150 LDA #result MOD 256 


160 STA 837 

170 LDA Zresult DIV 256 
180 STA 838 

190 LDA #1 

200 STA 839 

210 JSR iout 

220 RTS 

230 1 


240 NEXT pass% 
250 CALL start 
260 PRINT !result 
270 END 


>RUN 
1234 


EQUD 1234 
ғаш Ü 


\ 
\ 


a е ыны a GPO МӨ” Шым uut ut 


set up value 
result 


zeroise loop counter 


next byte of constant 
set up IWA 

bump Loop counter 
end of Loop 

no - back 

set up 

857 

set ир 

858 

set ир 

939 

call iout 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: iplus 

: IWA = МА + integer variable 
BASIC 1 address : 
: &9C5B 


&9C36 


: IWA contains a 32 bit integer 

: &04,&05 lo,hi point to integer variable 
: X = #&4 

: IWA added 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 50 microseconds 
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iplus demonstration (BASIC 2) 


0 iplus = 89С5В 
10 DIM mc% 200 


20 FOR pass% = 0 TO 2 STEP 2 


30 PZ = тс# 
40 СОРТ pass% 


50 .constant EQUD 26 
60 .adder EQUD 3 
70 .result EQUD 0 
80 .start 

90 LDX #0 

100 .Loop1 

110  LDA constant,X 

120 STA &2A,X 

130 INX 

140  CPX #4 

150 ВСС Loop1 

160 LDA Zadder MOD 256 
170 STA &4 

180 LDA #adder DIV 256 
190 STA &5 

200  LDX #4 

210  JSR iplus 

220 LDX #0 

230 .Loop2 

240 LDA &2A,X 

250 STA result,X 

260 INX 

270  CPX #4 

280 ВСС Loop2 

290 RTS 

300 1 


310 NEXT pass% 
320 CALL start 
330 PRINT result 
340 END 


>RUN 
29 


жет 


жет 


ИНЕ Е Р Же м, ilo Mu э". 


we gt. E со 


set up value 
result 
zeroise loop counter 


get next byte of constant 
save in IWA 

bump Loop counter 

end of Loop ? 

no - back 

get Lo address of adder 
save in &4 

get hi address of adder 
save in &5 

set X 

call iplus 

zeroise loop counter 


get next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 
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subroutine name : ipos 

function : Make IWA positive 

BASIC 1 address : &AD94 

BASIC 2 address : &AD71 

entry conditions : IWA contains a 32 bit integer 


comments : Ifthe IWA contains a negative integer, it is 
complemented. Else it is unchanged. 


exit status : IWA always positive 
: А destroyed 
: X unchanged 
: Ү destroyed 
: P destroyed 


typical timing : 17 or 27 microseconds (17 if already +ve) 
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ipos demonstration (BASIC 2) 


0 ipos = &AD71 
10 DIM mc% 200 


20 FOR pass% = 0 TO 2 STEP 2 


30 PZ = тс# 
40 СОРТ pass% 


50 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 


240 1 


LDX #0 


. Loop1 


LDA constant,X 
STA &2A,X 

INX 

CPX #4 

BCC Loop1 

JSR ipos 

LDX #0 


. Loop2 


LDA &2A,X 
STA result,X 
INX 

CPX #4 

BCC Loop2 
RTS 


250 NEXT pass% 
260 CALL start 
270 PRINT !result 
280 END 


>RUN 


1234 


.constant EQUD -1234 \ 
.result EQUD O 
-start 


\ 


set up value 
result 


zeroise loop counter 


get next byte of constant 
save in IWA 

bump Loop counter 

end of Loop ? 

no - back 

call ipos 

zeroise loop counter 


get next byte of IWA 
save in result 

bump Loop counter 
end of Loop ? 

no - back 
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subroutine name 


function 


BASIC 2 address 
entry conditions 


comments 


exit status 


typical timing 


: ismall 

: IWA = 256*Ү + A 
BASIC 1 address : 
: «АЕЕА 


&AF19 


: Ү and А set up appropriately 
: This is a handy way ofinitialising the IWA with 


small numbers up to 16 bits in length 


: IWA set up 

: А destroyed 
: X unchanged 
: Ү destroyed 
: P destroyed 


: 20 microseconds 
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ismall demonstration (BASIC 2) 


0 ismall = ФАЕЕА 

10 DIM mc% 200 

20 FOR pass% = 0 TO 2 STEP 2 
30 P% = mc% 

40 СОРТ pass% 


50 .result EQUD 0 \ result 

60 .start 

70 LDY #813 \ set up constant = #1388 
80 LDA #888 \ = 5000 

90 JSR ismall \ ismall to put 256*Y + A in IWA 
100 LDX #0 XV zeroise Loop counter 

110 .Loop 

120 LDA &2A,X \ get next byte of IWA 
130 STA result,X \ save in result 

140 INX \ bump Loop counter 

150  CPX #4 \ end of Loop ? 

160 ВСС Loop \ no - back 

170  RTS 

180 1 


190 NEXT pass% 
200 CALL start 
210 PRINT !result 
220 END 


>RUN 
5000 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


Comments 


exit status 


typical timing 


: itest 

: Test integer variable = > or < IWA 

: &9A85 

: &OAAD 

: IWA contains a 32 bit integer 

: &04,&05 (lo,hi) point to integer variable 
: &27 must be #4 

: On exit, the status register is set so that 
: BEQ will work if variable 2 IWA 

: BCC will work if variable « IWA 

: BCS will work if variable > ог = IWA 


: These tests must be performed immediately, or 


alternatively PHP can be used to stack the status 
register for testing later in the program 


: IWA destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P see comments above 
: &4,5 stepped by 4 


: 58 microseconds 
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itest demonstration (BASIC 2) 


10 itest = &9AAD : DIM mc% 200 
20 FOR pass% = 0 TO 2 STEP 2 
30 PX = mc% : СОРТ pass% 


40 .con1 EQUD 1234 \ set up value 1 
50 .con2 EQUD 2345 \ set up value 2 
60 .start LDX #0 XV zeroise [оор counter 
70 .loop1 LDA con2,X \ get next byte of con2 
80 STA &2A,X \ save in IWA 

90 INX \ bump Loop counter 
100 CPX #4 \ end of Loop ? 

110 BCC Loop1 \ no - back 

120 LDA #con1 MOD 256 \ get LSB of con1 
130 STA &4 \ save іп &4 

140 LDA #con1 DIV 256 \ get MSB of con1 
150 STA &5 \ save іп 85 

160 LDX #4 \ set for 

170 STX &27 \ integer variable type 
180 JSR itest \ call itest 

190 BEQ equals \ if = 

200 BCS more \ if > 

210 LDA #3 \ if < 

220 JMP set70 \ set &70 

230 .equals LDA #1 \ if = 

240 JMP set70 \ set &70 

250 .more LDA #2 \ if > 

260 .set70 STA &70 N &70 set 

270 RTS:1 

280 NEXT pass% : DIM test$(3) 

290 test$(1) = " is equal to " 

300 test$(2) - " is greater than " 

310 test$(3) = " is Less than " 

320 CALL start 

330 PRINT !con1;test$(?&70) ;!con2 

340 END 

>RUN 


1234 15 Less than 2345 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: izero 

: Make IWA zero 
BASIC 1 address : 
: &AECA 


&AEF9 


: none 
:IWA = 0 

: А destroyed 
: X unchanged 
: Y destroyed 
: P unchanged 


: 25 microseconds 
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ігего demonstration (BASIC 2) 


0 
10 
20 
50 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 


>RUN 


ігего = ФАЕСА 
DIM mc% 200 
FOR passX - 0 TO 2 STEP 2 


PZ = mc% 
COPT pass% 
-result EQUD 0 \ result 
-Start 
JSR ігего \ call izero 
LDX #0 \ zeroise Loop counter 
. Loop 
LDA &2A,X \ get next byte of IWA 
STA result,X \ save in result 
INX \ bump Loop counter 
CPX #4 \ end of loop ? 
BCC Loop \ no — back 
RTS 
1 
NEXT pass% 
CALL start 
PRINT !result 
END 
0 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 
Comments 


exit status 


typical timing 


: izpin 

: Copy integer variable in zero page to IWA 
: &AF85 

: &AF56 

: X points to zero page integer variable 
: &00,X to &03,X copied into IWA 

: IWA set up 

: А destroyed 

: X unchanged 

: Y unchanged 

: P destroyed 


: 27 microseconds 
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izpin demonstration (BASIC 2) 


0 izpin = &AF56 

10 DIM mc% 200 

20 FOR pass% = 0 TO 2 STEP 2 
30 PZ = тс# 

40 СОРТ pass% 


50 .result EQUD 0 \ result 

60 .start 

70 LDX #&70 \ point X to &70 

75 \ (contains 5000) 

80  JSR izpin \ call izpin 

90  LDX #0 XV zeroise Loop counter 
100 .(оор 

110 LDA &2A,X \ get next byte of IWA 
120 STA result,X \ save in result 

130 INX \ bump Loop counter 
140  CPX #4 \ end of Loop ? 

150  BCC Loop \ no - back 

160  RTS 

170 1 


180 NEXT pass% 
190 1870 = 5000 
200 CALL start 
210 PRINT !result 
220 END 


>RUN 
5000 
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subroutine name 


function 


BASIC 1 address : 
: &BE44 


BASIC 2 address 
entry conditions 
comments 


exit status 


typical timing 


: izpout 


: Copy IWA to integer variable in zero page 


&BE5C 


: X points to zero page integer variable 
: IWA copied to &00,X to &03,X 

: IWA unchanged 

: А destroyed 

: X unchanged 

: Y unchanged 

: P destroyed 

: &00,X to &03,X set up 


: 26 microseconds 
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izpout demonstration (BASIC 2) 


0 izpout = &BE44 

10 DIM mc% 200 

20 FOR pass% = O TO 2 STEP 2 
30 P% = mc% 

40 [COPT pass% 

50 .constant EQUD 5000 
60 .start 

70 LDX #0 

80 .Loop 


90  LDA constant,X 
100 STA %2А,Х 

110 INX 

120 CPX #4 

130 BCC Loop 

140  LDX #&70 

150 JSR izpout 

160  RTS 

170 1 


180 NEXT pass% 
190 CALL start 
210 PRINT !&70 
220 END 


>RUN 
5000 


\ 


\ 


ғ не” өтт PE gh 


set up value 
zeroise loop counter 


get next byte of constant 
save in IWA 

bump Loop counter 

end of Loop ? 

no - back 

point izpout to 870 

call izpout 
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3 FLOATING POINT 
NUMBERS 


Although floating point format can be used to store integers, 
especially large integers outside of the range of 32-bit integer 
format, it is essentially designed to handle fractions. 


3.1 Floating Point Variables 


The BBC BASIC interpreter recognises a floating point variable by 
the absence of either a ‘%’ or a ‘$’ at the end of its name. Each 
floating point variable occupies five bytes (40 bits). The number 
itself is held in the last four bytes (called the mantissa). The first 
byte (called the exponent) defines the position of the binary point. 
In other words, it defines the end of the integral part of the number 
and the start of the fractional part. 


Consider the decimal number 7.125. The binary equivalent of this 
number is 0111.0010 = (0*8) + (1*4) + (1*2) + (1*1) + (0*1/2) + 
(0*1/4) + (1*1/8) + (0*1/16). To obtain its floating point 
representation, the mantissa is simply written down as a string of 
32 bits, aligned at the most significant end, with the binary point 
omitted. 


mantissa 


0111 0010 ОООО 0000 0000 0000 0000 0000 
872 &00 &00 &00 


Had the binary point been included, it would have been positioned 
after the fourth bit in the mantissa. Thus the exponent is 4 in this 
case. 


exponent mantissa 


0000 0100 0111 0010 0000 0000 0000 0000 0000 0000 
&4 &72 &00 &00 &00 


It will be seen that the exponent and mantissa above are not the 
only ones that represent the number 7.125. Consider the following: 


exponent mantissa 


0000 0101 0011 1001 0000 0000 0000 0000 0000 0000 


&5 &39 &00 &00 &00 
0000 0110 0001 1100 1000 0000 0000 0000 0000 0000 
&6 &1C &80 &00 &00 
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When a zero is shifted into the most significant bit of the mantissa, 
the mantissa is effectively divided by 2. To compensate for this, the 
exponent is incremented by one, effectively multiplying the 
number by 2. The decimal analogy is that 7.125 could equally well 
be expressed as 71.25 tenths or 712.5 hundredths. Clearly it would 
be difficult to work with floating point numbers if each number 
could have во many different representations. 


Fortunately there is a rule which standardises the format of floating 
point numbers. The rule is called ‘Normalisation’. In a normalised 
floating point number, the most significant bit of the mantissa is 
always a 1. This is achieved by shifting the mantissa left, bit-by-bit, 
until all leading zeroes have disappeared. Since each leftward shift 
multiplies the mantissa by 2, the exponent must be reduced by one 
each time. Thus in its normalised form, the floating point 
representation of 7.125 is: 


exponent mantissa 


0000 0011 1110 0100 0000 0000 0000 0000 0000 0000 
&3 &E4 &00 &00 &00 


BBC BASIC expects all floating point numbers to be normalised. It 
makes use of this fact to handle the sign of a floating point number. 
If the sign of the number is positive, it changes the most significant 
bit of the mantissa to 0. This is simply a ruse to avoid holding the 
sign separately and hence to minimise the amount of memory 
needed to store a floating point number. In fact the number above 
represents —7.125 and +7.125 is: 


exponent mantissa 


0000 0011 0110 0100 0000 0000 0000 0000 0000 0000 
&3 &64 &00 &00 &00 


BBC BASIC adds &80 to the exponent. This is purely a device to 
assist in processing floating point numbers. Thus in BBC BASIC, a 
floating point variable set to +7.125 actually contains: 


exponent mantissa 


10000011 0110 0100 0000 0000 0000 0000 0000 0000 
885 864 &00 &00 &00 


It will be seen that, to convert a positive number to negative, it is 
only necessary to add &80 to the most significant byte of the 
mantissa. 
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Thus there are four stages in the process of converting a number to 
floating point format: 


а) Write down the mantissa in binary and set the exponent to the 
value which fixes the position of the implied binary point. 


b) Normalise the number, ensuring that the exponent is adjusted 
appropriately for each bit shifted. 


c) Ifthe sign of the original number was positive, change the most 
significant bit of the mantissa to zero. 


d) Add &80 to the exponent. 


It remains now to demonstrate that this process works for purely 
fractional numbers. Consider the decimal number -0.375. As this 
represents -3/8 which is the same as -(1/4 + 1/8), the binary 
equivalent of this number is -.0110. It is easy to write down the 
mantissa in its unnormalised binary form: 


mantissa 


0110 0000 ОООО 0000 0000 0000 0000 0000 
860 &00 &00 &00 


The exponent of this unnormalised number is zero, since the 
implied binary point comes immediately in front of the most 
significant bit of the mantissa. To normalise the mantissa, it must 
be shifted to the left until all leading zeroes have been removed, 
reducing the exponent by one each time. Consequently, the 
exponent is –1, to which must be added &80 as before. In 
hexadecimal terms: 


-0.375 = &7F &CO &00 &00 &00 


One last contrivance is employed in BBC BASIC. The floating point 
representation of 0.0 does not follow the rules. It is simply stored as 
five bytes of zeroes. 


Some floating point numbers are tabulated below: 


0 = &00 &00 &00 &00 %00 
*1 - &81 &00 &00 &00 &00 -1 - &81 &80 &00 &00 &00 
+2 = 882 &00 &00 &00 &00 -2 - &82 &80 &00 &00 &00 
*5 - &82 &40 &00 &00 &00 -3 = 882 &CO &00 &00 %00 
*4 - 885 &00 &00 &00 &00 -4 = 885 880 &00 &00 &00 
%5 - 885 820 &00 &00 &00 -5 = 885 &AO 800 800 &00 
+6 = 885 840 800 800 &00 -6 = 885 &CO 800 800 %00 
+7 = 885 860 800 800 &00 -7 = 885 8Е0 800 800 &00 
+8 = 884 800 &00 800 %00 -8 = 884 880 &00 800 %00 
+9 = 884 810 800 800 &00 -9 = 884 890 800 800 %00 
+10 = 884 820 &00 500 8000 -10 = 884 «А0 200 800 800 
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3.2 Integer versus Floating Point 


Integer numbers have two big advantages over their floating point 
counterparts: accuracy and speed. 


In all the floating point examples so far presented in this book, the 
fractional part of the number has always translated into an exact 
number of halves, quarters, eighths, sixteenths and so on. Numbers 
like this are referred to as ‘machine numbers’ because they can be 
held exactly within a computer. By no means can this be said of all 
numbers. For example, the fraction 1/5 cannot be held exactly. 
However many fractional bits are included, the resulting floating 
point number remains an approximation, albeit a good one, to the 
original fraction. Numbers that have a large integral part have less 
bits available for the fractional part and tend to suffer more from 
fractional inaccuracy. An example of floating point inaccuracy is 
shown in the following BASIC code: 


10 J-0 

20 FOR I% = 1 TO 40 
30 J = J + 0.2 

40 NEXT 

50 PRINT J 

60 END 


>RUN 
7.99999999 


Integer arithmetic is not only accurate, it is faster than floating 
point. This simply means that it is less complicated and requires 
less code in the BASIC interpreter. 


However, integer numbers have two disadvantages compared to 
floating point numbers. The first is obvious; they cannot be used to 
represent fractions. The second is that they cannot handle such a 
wide range of numbers as floating point (from 1.7*10 to the power 
-39 up to 1.7*10 to the power 38). 


Overall the advantages of using integers are so great, that providing 
the numbers to be used fall within the range of integer numbers, 
they should be used if at all possible. This is especially true for 
financial applications, where loss of accuracy in floating point can 
render a program useless. Except for very large sums of money, 
financial data should be held in pence. There is then no fraction to 
consider (assuming the demise of the halfpenny). The decimal 
point can always be inserted into an ASCII field when printing 
reports. 
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This is an example of a technique known as scaling. Because 
financial data has only two digits following the decimal point, all 
numbers are scaled up by a factor of 100 to remove the fractional 
part altogether. The technique can be applied to other situations. 


Floating point comes into its own in mathematical programs where 
total accuracy is neither expected nor required, such as when 
drawing curves. 


3.3 Floating Point Work Areas 


BBC BASIC does all its floating point arithmetic in two working 

areas in zero page. Ав with the IWA, these memory areas are not 
dedicated to floating point and may be re-used for quite different 
purposes. 


Each of the working areas is eight bytes long (not five as might have 
been expected). The extra 3 bytes are for the following purposes: 


а) a sign byte. In the five byte variable format, the sign is 
contrived by zeroising the most significant bit of the mantissa 
for positive numbers. In eight byte format, the most significant 
byte of the mantissa is copied into a sign byte, and the most 
significant bit of the mantissa is restored to 1. 


b) a rounding byte. An extra byte is tacked onto the end to extend 
the precision of arithmetic. This extra byte is used to round the 
preceding mantissa when required. 


c) anoverflow byte. This byte exists to trap errors, such as those 
which might occur when the result of a multiplication is a 
number too big to handle. 


The five byte format unpacks into the eight byte format as follows: 


5 byte 8 byte 

exponent byte 0 —— > byte 2 exponent 
mantissa-1 byte 1 ------ > byte 0 sign 
mantissa-1 byte 1 OR 8880 byte 3 mantissa-1 
mantissa-2 byte 2 ---- > byte 4 mantissa-2 
mantissa-3 byte 5 ------ > byte 5 mantissa-3 
mantissa-4 byte 4 ------ > byte 6 mantissa-4 

zero byte 1 overflow 

zero byte 7 rounding 


59 


For example, +1.0 in the two formats is: 


5 byte 8 byte 
exponent 881 800 sign 
mantissa-1 800 800 overflow 
mantissa-2 800 881 exponent 
mantissa-3 800 880 mantissa-1 
mantissa-4 800 800 mantissa-2 


800 mantissa-3 
800 mantissa-4 
800 rounding 


The two floating point areas used by BASIC are &2E % &35, and 
&3B to &42 inclusive. These will be referred to as Floating Point 
Work Area A or FWA, and Floating Point Work Area B or FWB, 
respectively. They consist of: 


sign &2E &3B 
overflow &2F &3C 
exponent 850 &3D 
mantissa-1 &31 &3E 
mantissa-2 &32 &3F 
mantissa-3 855 &40 
mantissa-4 854 &41 
rounding 855 &42 
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3.4 Defining Floating Point Constants 


Neither BASIC 1 nor BASIC 2 has a directive to allow definition of 
a floating point constant to the assembler. Once again it is 
necessary to invent a pseudo-directive. This one is called EQUF. It 
relies on two facts. Firstly it uses a BASIC variable, Z, and the 
address of the look-up table for variables starting with the letter ‘Z’ 
is held іп %4В4,%4В5 (lo,hi). Secondly, the pseudo-directive 
assumes that the actual data in Z will be found 3 bytes into this 
look-up table. This will be true so long as Z is the first variable 
which has an entry in this look-up table. To ensure this, make 
certain that no other BASIC variables start with the letter 'Z'. This 
caution should be unnecessary, because the techniques advocated 
discourage hybrid assembly language/BASIC programs. 


10 DIM mc% 100 : FOR pass% = 0 TO 2 STEP 2 
20 PX = mc% : СОРТ pass% 

30 .constant 

40 ОРТ FNEQUF(1.0) 

50 1 

60 NEXT pass% 

TO END 

80 DEF FNEQUF(Z) 

90 1% = 3 + ?&4B4 + 256*?84B5 

100 FOR JZ = 1 TO 5 


110 ЭРХ = 212 
120 P% = РХ + 1 
130 IX = 14 +1 
140 МЕХТ 


150 = pass% 


A particularly useful feature of these pseudo-directives, is that they 
can be used to evaluate expressions, providing the terms are also 
literals. For example: 


30 .constant 
40 ОРТ FNEQUF(3XxSIN(PI/4)) 
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3.5 Floating Point Routines Summary 


The following tables summarise floating point routines available 
within the BASIC ROM. Conversion routines, e.g. floating point to 
ASCII, are the subject of a later chapter. So too are routines 
handling trigonometric functions, square roots, logarithms etc. 


Name BASIC 1 BASIC 2 Function 
address address 

aclear  &A691 84686 FWA - 0 

acomp &AD99  &AD7E FWA - -ҒМА 

acopyb &A4E4 — &AADC FWA - FWB 

adiv &A68B &АбАР” FWA - fp var / ҒМА normalised 
and rounded 

adiv10  &A23E &A24D FWA = FWA / 10 unnormalised 
and unrounded 

aminus  &A50B  &A4FD FWA = fp var — FWA normalised 
and rounded 

amult &A661 &A656 FWA - ҒМА ж fp var normalised 
and rounded 

amult1 | &A611 &А606 FWA = FWA * fp var 
unnormalised and unrounded 

amult10 &A1E5  &A1F4 FWA = FWA ж 10 unnormalised 
and unrounded 

anorm &A2F4  &A303 normalise FWA 

aone &AGA4 &A699 FWA = 1 

apack &A37E  &A38D pack FWA into fp var 

apack1 &A376 &A385 pack FWA into &46С to &470 

apack2 &A36E &A37D pack FWA into 8471 to &475 

apack3 &A372  &A381 pack FWA into 8476 to &47A 

aplus &A5OE  &A500 FWA = FWA + fp var normalised 
and rounded 

aplusb &A513 &А505 FWA = FWA + FWB normalised 
and rounded 

aplus1 ----- &A50B FWA = FWA + FWB normalised 
and unrounded 

arecip  &A6BO &A6A5 FWA = 1 / FWA normalised 
and rounded 

around %А667 &A65C round FWA 

asign &A1CB — &ATDA get sign of FWA 

aswap &A4DE  &A4D6 swap fp var and FWA 

atest &9A37 %9А5Ғ test fp var against FWA 

aunp &A3A6 &A3B5 unpack fp var into FWA 

aunp1 &A3A3 &A3B2 unpack &46C to &470 into FWA 

bclear 84465 &A453 FWB = 0 

bcopya  &A20F  &A21E FWB = FWA 

bunp &A33F &A34E unpack fp var into FWB 
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3.6 Floating Point Routines Description 


The following pages show how to use each of the floating point 
routines in the table. Set up procedures are very simple. Routines 
involving only FWA and/or FWB need no setting up, other than to 
load FWA/FWB prior to the call. Routines which reference а 
floating point variable should point &4B,&4C (lo,hi) to that variable. 


For this reason, demonstration programs are not supplied for each 
routine. In place of these, there is a single program which can be 
incorporated into a user program, which: 


a) provides a standard interface for all floating point arithmetic. 
b) can be assembled under BASIC 1 or BASIC 2. 


c) the machine code derived will run under either BASIC 1 or 
BASIC 2. 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: aclear 

: FWA = zero 
BASIC 1 address : 
: &A686 


&A691 


: none 
: FWA = zero 

: FWB unchanged 
: А destroyed 

: X unchanged 

: Ү unchanged 

: P destroyed 


: 25 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: acomp 
: FWA = -FWA 
BASIC 1 address : 
: &AD7E 


&AD99 


: none 
: FWA = -FWA 

: FWB unchanged 
: A destroyed 

: X unchanged 

: Y unchanged 

: P destroyed 


: 34 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: acopyb 

: FWA = FWB 
BASIC 1 address : 
: &A4DC 


&AAEA 


: none 
: FWA = FWB 

: FWB unchanged 
: А destroyed 

: X unchanged 

: Y unchanged 

: P destroyed 


: 36 microseconds 
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subroutine name : айу 
function : FWA = fp var / FWA normalised, rounded 
BASIC 1 address : &A68B 
BASIC 2 address : A6AD 
entry conditions : &4B,&4C (lo,hi) point to fp var 
exit status : FWA = quotient 
: FWB destroyed 
: A destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 
error reports : division by zero 


typical timing : 1545 microseconds 
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subroutine name : adiv10 
function : FWA = FWA / 10 normalised, rounded 
BASIC 1 address : &A23E 
BASIC 2 address : &A24D 
entry conditions : none 
exit status : FWA = FWA / 10 
: FWB destroyed 
: А destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 


typical timing : 360 microseconds 
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subroutine name : aminus 
function : FWA = fp var - FWA normalised and rounded 
BASIC 1 address : &A50B 
BASIC 2 address : &A4FD 
entry conditions : &4B,&4C (lo,hi) point to fp var 
exit status : FWA = result 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


typical timing : 254 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: &A656 


BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: amult 


: FWA = fp var * FWA normalised and rounded 


&A661 


: &4B,&4C (lo,hi) point to fp var 
: FWA = result 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: too big 


: 1581 microseconds 
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subroutine name : 
: FWA = fp var * FWA unnormalised and 


function 


BASIC 1 address : 
BASIC 2 address : 


entry conditions 


exit status 


error reports 


typical timing 


amult1 


unrounded 
&A611 
&A606 


: &4B,&4C (lo,hi) point to fp var 
: FWA = result 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: too big 


: 1508 microseconds 
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subroutine name 


function 


BASIC 1 address : 
BASIC 2 address : 


entry conditions 


exit status 


typical timing 


: amult10 
: FWA = 10 * FWA unnormalised and 


unrounded 
&A1E5 
&A1F4 


: none 
: FWA = result 

: FWB destroyed 
: А destroyed 

: X unchanged 

: Ү unchanged 

: P destroyed 


: 171 microseconds 


72 


subroutine name 


function 


BASIC 1 address : 
: & A303 


BASIC 2 address 
entry conditions 


exit status 


typical timing 


: anorm 


: FWA = FWA normalised 


&A2F4 


: none 
: FWA = result 

: FWB unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 27 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: aone 
: FWA=1 
BASIC 1 address : 
: &A699 


&АбА4 


: попе 
: FWA = 1 

: FWB unchanged 
: А destroyed 

: X unchanged 

: Y destroyed 

: P destroyed 


: 37 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: apack 

: fp var = FWA 
BASIC 1 address : 
: &A38D 
: &4B,&4C (lo,hi) point to fp var 
: FWA unchanged 
: FWB unchanged 
: A destroyed 


&A37E 


: X unchanged 
: Y destroyed 
: P destroyed 


: 46 microseconds 


75 


subroutine name : apack1 
function : &46C to &470 = FWA 
BASIC 1 address : &A376 
BASIC 2 address : &A385 
entry conditions : none 
exit status : FWA unchanged 
: FWB unchanged 
: A destroyed 
: X unchanged 
: Y destroyed 
: P destroyed 
: &4B,&4C destroyed 


typical timing : 51 microseconds 
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subroutine name : apack2 
function : &471 to &475 = FWA 
BASIC 1 address : &A36E 
BASIC 2 address : &A37D 
entry conditions : none 
exit status : FWA unchanged 
: FWB unchanged 
: А destroyed 
: X unchanged 
: Y destroyed 
: P destroyed 
: &4B,&4C destroyed 


typical timing : 53 microseconds 
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subroutine name : apack3 
function : &476 to &A7A = FWA 
BASIC 1 address : &A372 
BASIC 2 address : &A381 
entry conditions : none 
exit status : FWA unchanged 
: FWB unchanged 
: А destroyed 
: X unchanged 
: Y destroyed 
: P destroyed 
: &4B,&4C destroyed 


typical timing : 53 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: &A500 


BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: aplus 


: FWA = fp var + FWA normalised and rounded 


&A50E 


: &4B,&4C (lo,hi) point to fp var 
: FWA = result 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: too big 


: 246 microseconds 
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subroutine name 


function 


BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: aplusb 

: FWA = FWA + FWB normalised and rounded 
BASIC 1 address : 
: &A505 


&A513 


: none 
: FWA = result 

: FWB destroyed 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: too big 


: 58 microseconds 
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subroutine name 


function 


BASIC 1 address : 
BASIC 2 address : 


entry conditions 


exit status 


error reports 


typical timing 


: aplus1 
: FWA = FWA + FWB normalised and 


unrounded 


: none 
: FWA = result 

: FWB destroyed 
: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: too big 


: 41 microseconds 
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subroutine name : arecip 
function : FWA = 1 / FWA normalised and rounded 
BASIC 1 address : «А6В0 
BASIC 2 address : &A6A5 
entry conditions : none 
exit status : FWA = result 
: FWB destroyed 
: A destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 
: &476 to &47A destroyed 


typical timing : 1619 microseconds 


82 


subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: around 
: FWA = FWA rounded 
BASIC 1 address : 
: &A65C 


&A667 


: none 
: FWA = result 

: FWB unchanged 
: А destroyed 

: X unchanged 

: Ү unchanged 

: P destroyed 


: 22 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: «АЛПА 


BASIC 2 address 
entry conditions 


comments 


exit status 


typical timing 


: asign 


: A register denotes sign of FWA 


&A1CB 


: none 
: A register set as follows: 

: = 0 if FWA is zero 

: = 1 if FWA is +ve 

: = -ve if FWA is -ve otherwise 
: FWA unchanged 

: FWB unchanged 

: A see above 

: X unchanged 

: Y unchanged 

: P destroyed 


: 24 microseconds 
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subroutine name : aswap 
function : FWA = fp var and fp var = FWA 
BASIC 1 address : &A4DE 
BASIC 2 address : &A4D6 
entry conditions : &4B,&4C (lo,hi) point to fp var 
exit status : FWA = fp var 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


typical timing : 121 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: &9A5F 


BASIC 2 address 
entry conditions 


comments 


exit status 


typical timing 


: atest 


: test fp var against FWA 


&9A37 


: &4B,&4C (lo,hi) point to fp var 

: on exit, the P register is set up such that: 
: BEQ works if fp var = FWA 

: BCC works if fp var « FWA 

: BCS works if fp var > or = FWA 
: FWA destroyed 

: FWB destroyed 

: A destroyed 

: X destroyed 

: Y destroyed 

: P see above 


: 78 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: aunp 
: FWA = fp var 
BASIC 1 address : 
: &A3B5 
: &4B,&4C (lo,hi) point to fp var 
: FWA = fp var 

: FWB unchanged 
: A destroyed 


&A3A6 


: X unchanged 
: Y destroyed 
: P destroyed 


: 49 microseconds 
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subroutine name 


function 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: aunp1 
: FWA = &46C to &470 
BASIC 1 address : 
: &A3B2 


&A3A3 


: none 
: FWA = result 

: FWB unchanged 
: А destroyed 

: X unchanged 

: Y destroyed 

: P destroyed 


: 62 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


: bclear 

: FWB=0 

: &A463 

: &A453 

: none 

: FWA unchanged 
: FWB=0 

: А destroyed 
: X unchanged 
: Y unchanged 
: P destroyed 


: 25 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


: bcopya 

: FWB = FWA 
: &A20F 

: &A21E 

: none 

: FWA unchanged 
: FWB = FWA 
: А destroyed 
: X unchanged 
: Ү unchanged 
: P destroyed 


: 36 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


: bunp 

: FWB = fp var 
: &A33F 

: &АЗ4Е 

: &4B,&4C (lo,hi) point to fp var 
: FWA unchanged 
: FWB = result 

: A destroyed 

: X destroyed 

: Y unchanged 

: P destroyed 


: 51 microseconds 
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3.7 Floating Point Interface Program 


This program will perform addition, subtraction, multiplication or 
division in floating point. It is intended that this code is 
incorporated into a user program. It may be assembled under either 
BASIC 1 or BASIC 2. When executed it will detect which BASIC is 
present and operate accordingly. 


On entry to floatsub, the following fields must be set up: 


870,871 Clo,hi) point to argument 1 
&72,8&73 (lo,hi) point to argument 2 
878,875 (lo,hi) point to argument 3 
&76 = function required 


if &76 = 0 argument 3 = argument 1 + argument 2 
if &76 = 1 argument 3 = argument 1 - argument 2 
if &76 = 2 argument 3 = argument 1 * argument 2 
if &76 = 3 argument 3 = argument 1 / argument 2 


АП mathematical functions return results into argument 3. АП 
results are fully normalised and rounded. Signs are handled 
correctly throughout. All registers are returned uncorrupted. 


Standard Floating Point Interface 


10 DIM mcZ 500 

20 FOR pass% = 0 TO 2 STEP 2 
30 PZ = тс# 

40 СОРТ pass% 

50 .basic2 \ BASIC 2 
60 .aplus JMP &A500 

70 .aminus JMP &A4FD 

80 .amult JMP &A656 

90 .adiv JMP &AÓAD 

100 .араск ./МР &A38D 

110 .aunp JMP &A3B5 


120 .basic1 \ BASIC 1 

130 JMP &A50E 

140 JMP &A50B 

150 JMP &A661 

160 JMP &A68B 

170 JMP &A37E 

180 JMP &A3A6 

190 .floatsub 

200 PHP \ save P reg 

210 PHA \ save А reg 

220 TXA \ save X reg 

230 PHA 

240 TYA \ save Y reg 

250 PHA 

260 LDA &8015 \ get BASIC year 
270 CMP #&32 \ is it 1982 ? 
280 BEQ skipmove \ yes - skipmove 
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290 LDX #0 
500 .move 


et 


zeroise loop counter 


310 LDA basic1,X \ next byte of BASIC 1 addresses 
320 STA basic2,X \ overwrite BASIC 2 counterpart 
330 INX \ bump loop counter 

340 CPX #18 \ end of loop ? 

350 BCC move \ no - move 

360 .skipmove 

370 LDA &72 \ get LSB of argument 2 
380 STA &4B \ save in &4B 

390 LDA &73 \ get MSB of argument 2 
400 STA &4C \ save in &4C 

410 JSR aunp \ unpack argument 2 into FWA 
420 LDA &70 \ get LSB of argument 1 
430 STA &4B \ save in &4B 

440 LDA &71 \ get MSB of argument 1 
450 STA &4C \ save in &4C 

460 LDA &76 \ get function 

470 BEQ add \ if add 

480 CMP #1 \ test subtract 

490 BEQ subtract \ if subtract 

500 СМР #2 \ test multiply 

510 BEQ multiply \ if multiply 

520 CMP #3 \ test divide 

530 BEQ divide \ if divide 

540 BRK \ invalid 

550 .add 

560 JSR aplus \ do add 

570 JMP result \ output result 

580 .subtract 

590 JSR aminus \ do subtract 

600 JMP result \ output result 

610 .multiply 

620 JSR amult \ do multiply 

630 JMP result \ output result 

640 .divide 

650 JSR adiv \ do divide 

660 .result 

670 LDA &74 \ get LSB of result 

680 STA &4B \ save in &4B 

690 LDA &75 \ get MSB of result 

700 STA &4C \ save in &4C 

710 JSR apack \ output result 

720 .restore 

730 PLA \ restore 

740 TAY N Y reg. 

750 PLA \ restore 

760 ТАХ \ X reg. 

770 PLA \ restore А reg 

780 PLP \ restore P reg 

790 RTS:1 


800 NEXT pass%:STOP 
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3.8 Floating Point Interface Program Tested 


The previous program can be tested by changing the BASIC 
statements. In the changes below, data values are entered into A 
and B (lines 900 and 910) and the program performs all the 
calculations. 


800 REM test all functions 
810 NEXT pass% 

820 DIM func$(4) 

830 func$(1) 
840 func$(2) 
850 func$(3) 
860 func$(4) 
870 ?&76 = -1 
880 REPEAT 
890 2876 = 2976 + 1 

900 А = -20.5 

910 B = -7.258 

920 C = 0 

930 I4 = 3 + 28482 + 256*?&&83 
940 2870 = IX MOD 256 


Hl H HI 
| + 


М + 


950 2871 = 1% DIV 256 
960 I% = 3 + 28484 + 256х28485 
970 2872 = 1% MOD 256 
980 2875 = 17 DIV 256 
990 IZ = 3 + 78486 + 256х28487 
1000 ?&74 = 1% MOD 256 


1010 2%75 = 17 DIV 256 

1020 CALL floatsub 

1030 PRINT A;func$(14?876);B;" = ";C 
1040 UNTIL 2876 - 5 


1050 END 

>RUN 
-20.3 + -7.258 = -27.558 
-20.3 - -7.258 = -13.042 
-20.3 ж -7.258 = 147.3374 
-20.3 / -7.258 = 2.79691375 
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4 CONVERSIONS 


Inevitably, assembly language programs have to deal with the 
problem of data conversion, particularly from binary to ASCII and 
vice-versa. The BASIC ROM contains a set of accessible conversion 
routines which are among the handiest of all the routines in the 
ROM. 


4.1 Conversion Work Areas 


The BASIC ROM does all its work with ASCII strings in an area of 
memory starting at &600. This will be called the String Work Area 
or SWA. There is another memory area, &36, which contains the 
length of the current string in the SWA. Another way of looking at 
this field is that the contents of &36, when added to #&600, point to 
the next available space in the SWA. 


The usual caution is necessary at this point. The area starting at 
&600 is not reserved for a dedicated purpose. It is also used as a 
variable parameter block. As for &36, its uses are legion. 


Another zero page location which plays an active role in some of 
these routines is &15. &15 controls the radix when converting a 
number to an ASCII string. If set to zero, the ASCII string represents 
a decimal number. If set to —1, it is hexadecimal. BASIC itself sets 
this field for the PRINT command. The latter setting is used if a 
tilde (~) appears іп the PRINT statement. 


During conversions to ASCII, parts of the print format field, 996, 
are important. Located from &400 to &403, this field controls print 
formatting as follows: 


&403 not applicable (STRS flag) 
&402 format number 


format 0 - general 
format 1 - exponential 
format 2 - fixed decimal 


&401 number of digits (exact meaning depends on format) 


format 0 = maximum number of digits which сап be 
printed before exponential format is used 
instead 

format 1 = number of digits + 1 that follows the ‘E’ 

format 2 = number of decimal places 


&400 width of print field 
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%400 does not affect string conversion, but rather is used by the 
BASIC PRINT command to work further on the converted string. 
&401 and &402 do affect string conversion. Format 0 is the default, 
representing typical BASIC formatting. Format 1 specifies that 
numbers are to be converted to exponential format e.g. 1000 would 
be converted to 1E3. Format 2 specifies that the number is to be 
converted to ASCII with a fixed number (to be found in &401) of 
decimal places. The maximum number of decimal places that can 
be specified is ten. If a number greater than 10 is placed in &401, 
the conversion routines default to ten. Similarly, if a format greater 
than 2 is specified, a default of zero is used. After BREAK, or at 
power up, BASIC initialises @% to format zero and a width of 10. 


4.2 Conversion Routines Summary 


Many of the routines below share common entry points. As the set 
up parameters are quite different in each case, it is easier to present 
them separately. Integer to floating point conversions (and vice- 
versa) are so straightforward that no further explanation is required. 
A demonstration program is supplied in the next section to show 
the use of ASCII conversion routines. 


Name BASIC 1 BASIC 2 Function 

address address 
ascnum &АСЅА  &AC34 ASCII to integer ог 

floating point 

fpascdec &9EDO &9ЕрЕ floating point to decimal ASCII 
fpaschex &9EDO &9ЕрЕ floating point to hex ASCII 
fpi1 &A3F2 %АЗЕ4 floating point to integer 1 
fpi2 &A40C — &ASFE floating point to integer 2 
jascdec &9EDO &9EDF integer to decimal ASCII 
jaschex &9EDO &9EDF integer to hex ASCII 
ifpa &A2AF  &A2BE integer to floating point 
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4.3 Conversion Routines Description 


subroutine name : ascnum 


function : the ASCII string in SWA is converted to either 
an integer placed in IWA or to a floating point 
number placed in FWA 


BASIC 1 address : &AC5A 

BASIC 2 address : «АСЗ4 

entry conditions : &36 contains length of string in SWA 
: SWA contains ASCII number 


comments : the routine places a binary zero at the end of 
SWA and steps &36 by 1 


: on exit A and &27 both reflect the type of 
conversion: 


: = #&40 for integer 

: = #&ЕЕ for floating point 
exit status : IWA = result (integer) 

: FWA = result (floating point) 

: FWB destroyed 

: SWA has zero appended 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &19,&1A,&1B destroyed 

: &45,&48,&49 destroyed 

: &27 see above 


typical timing : 1748 microseconds 
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subroutine name 


function 


BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


typical timing 


: fpascdec 


: the floating point number in FWA is converted 


to ASCII decimal and placed in SWA 


: &9EDO 

: &9EDF 

: &15 must be zero 

: Y must be #&FF 

: @% must be set as appropriate 


: the routine returns with &36 set to the length of 


the string 


: IWA destroyed 

: FWA destroyed 

: FWB destroyed 

: SWA = result 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &36 - length of string 


: 4878 microseconds 
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subroutine name 


function 


BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


typical timing 


: fpaschex 


: the floating point number in FWA is converted 


to ASCII hexadecimal and placed in SWA 


: &9EDO 

: &9EDF 

: &15 must be #&FF 

: Ү must be #&FF 

: @% must be set as appropriate 

: only the integer part of the number is converted 


: the routine returns with &36 set to the length of 


the string 


: IWA destroyed 

: FWA destroyed 

: FWB destroyed 

: SWA = result 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &36 - length of string 


: 582 microseconds 


99 


subroutine name : fpi1 


function : the floating point number in FWA is converted 
to integer and placed in IWA 


BASIC 1 address : &A3F2 
BASIC 2 address : &A3EA 
entry conditions : none 
comments : only the integer part of the number is converted 
exit status : IWA = result 
: FWA destroyed 
: FWB destroyed 
: SWA unchanged 
: А destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 


typical timing : 405 microseconds 
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subroutine name 


function 


BASIC 1 address 
BASIC 2 address 
entry conditions 


comments 


exit status 


typical timing 


: fpi2 


: the floating point number in FWA is converted 


to integer and placed in &31 to &34 


: &A40C 

: &АЗЕЕ 

: none 

: only the integer part of the number is converted 


: this routine сап be used instead of fpi1 when it 


is required to preserve IWA 


: IWA unchanged 

: FWA destroyed 

: FWB destroyed 

: SWA unchanged 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &31 to &34 = result 


: 387 microseconds 
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subroutine name : 


function 


BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


typical timing 


iascdec 


: the integer number in IWA is converted to 


ASCII decimal and placed in SWA 


: &9EDO 

: &9EDF 

: &15 must be zero 

: Y must be #&40 

: @% must be set as appropriate 


: the routine returns with &36 set to the length of 


the string 


: IWA destroyed 

: FWA destroyed 

: FWB destroyed 

: SWA = result 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &36 - length of string 


: 5369 microseconds 


102 


subroutine name : 


function 


BASIC 1 address 
BASIC 2 address 


entry conditions 


comments 


exit status 


typical timing 


iaschex 


: the integer number in IWA is converted to 


ASCII hexadecimal and placed in SWA 


: &9EDO 

: &9EDF 

: &15 must һе #&FF 

: Ү must be #&40 

: @% must be set as appropriate 


: the routine returns with &36 set to the length of 


the string 


: IWA destroyed 

: FWA destroyed 

: FWB destroyed 

: SWA = result 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: &36 - length of string 


: 216 microseconds 
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subroutine name 


function 


BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


: ifpa 


: the integer number in IWA is converted to 


floating point and placed in FWA 


: &A2AF 

: «А2ВЕ 

: none 

: IWA destroyed 
: FWA = result 

: FWB destroyed 
: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 149 microseconds 
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4.4 ASCII Conversion Demonstration 


The following program is written in BASIC 2. It performs the 
following tasks: 


a) 


b) 
с) 


а) 


е) 


It asks Ше user to supply а number. This сап be either integer 
or decimal, with or without a sign. 


It uses OSWORD 0 to read the user's number into the SWA. 


It converts this number to either binary integer or floating point 
using ‘ascnum’ 


It then converts the number back into fixed format ASCII 
decimal with 5 decimal places, overwriting the SWA. This 
result is then printed. 


note that the conversion process in d) references a routine 
called ‘numasc’ which represents either fpascdec or iascdec, 
this being controlled by the value in Y prior to the call. 


ascnum = &AC34:numasc = &9EDF 
osword = &FFF1:osbyte = &FFF4 
osnewl = &FFE7:oswrch = &FFEE 
DIM mc% 500 
FOR pass% = O TO 2 STEP 2 
PZ = mc% 
COPT pass% 
.msg EQUB 12 \ clear screen 
EQUB 31 \ cursor 
EQUB 1 \х=1 
EQUB 12 \у=12 
EQUS "enter your number >" \ message 
-msgl EQUB msgl - msg \ msg Length 
.osbuf f EQUW &600 \ point to SWA 
EQUB 20 \ max size 
EQUB &2A \ min value 
EQUB &39 \ max value 
-osbuffadd EQUW osbuff \ pointer to osbuff 
-Start 
LDX #0 \ zeroise Loop counter 
. Loopmsg 
LDA msg,X \ get next byte of msg 
JSR oswrch \ print it 
INX \ bump X 
CPX msgl \ end of msg ? 
BNE Loopmsg \ no — back 
\ get a number from keyboard 
.reply 
LDX osbuffadd \ X = to 
LDY osbuffadd+1 \Y=hi 
LDA #0 \ OSWORD 0 
JSR оѕмога \ to read reply 
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520 ВСС 
550 LDA 
340 JSR 
350 JMP 
360 .notesc 
370 STY 
380 JSR 
390 TAY 
400 LDA 
410 STA 
420 LDA 
430 STA 
440 LDA 
450 STA 
460 JSR 
470 JSR 
480 LDX 
490 .ploop 
500 LDA 
510 JSR 
520 INX 
530 CPX 
540 BCC 
550 JSR 
560 RTS 
570 1 


notesc 
#&ТЕ 
osbyte 
reply 


&36 
ascnum 


#0 

815 

#2 
8402 
#5 
8401 
numasc 
osnewl 
#0 


&600,X 
oswrch 


856 
ploop 
osnewl 


580 NEXT pass% 
590 CALL start 


600 END 
>RUN 


enter your number >1234.567 


1234.56 


700 


CC CCC лжештемжшжт 


CCC 
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not ESCAPE 
acknowledge 
ESCAPE 

try again 


set up reply Length 
convert to binary 
save variable type 
set decimal 

print 

set format 

=2 

set decimal places 
=5 

convert it back 
new Line 

zeroise loop counter 


get next byte of ASCII 
print it 

bump X 

end of print 

no - back 

new Line 


5 MATHEMATICAL 
FUNCTIONS 


The mathematical functions are surprisingly simple to use. No new 
ground has to be covered to explain their use. However, the way in 
which many of them work may be of some interest. It certainly has 
a bearing on the time they take to run. 


It might be expected that the BASIC ROM would contain tables of 
sines, cosines etc. This is not true. Tables would use up far too 
much memory. Most of the mathematical functions can be 
expressed as series. For example: 


а ажа ахаха ахахажа 
ехр(а)=1 + - + — + —— + —— Pan same dass etc. 
7 9% 32221  4x3x2x1 


This is an example of a convergent series. Each successive term in 
the series provides a value smaller than the previous term. Thus if 
enough terms are taken a good approximation results. The more 
terms that are taken, the longer it takes to execute; fewer terms 
reduce the accuracy of the final answer. In practice, therefore, the 
number of terms considered is a compromise between accuracy and 
execution time. 


In any event, these routines can never be fast to execute. It follows 
that if they are used in a loop with a large number of iterations, the 
effect on execution time is significant. The demonstration program 
in this chapter, which is only intended to show how to use the 
mathematics routines, is an example of slow circle drawing. The 
slowness is due to the fact that both sine and cosine functions are 
used within a loop. 
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5.1 Mathematical Functions Routines Summary 


АП of the functions summarised below, bar pi, work in the same 
way. The floating point number on which the function is to operate 
is placed in FWA. After the routine has been executed, the required 
result is to be found in FWA. Note that acs requires two subroutine 
calls, one immediately after the other. Note also that each function 
is equivalent to the BASIC function with the same name (but in 
upper case letters). 


Name BASIC 1 BASIC 2 Function 
address address 
acs 1) &A8CF &A8DD FWA = acs (FWA) 
2) &A929 &A927 


asn &A8CF &A8DD FWA = asn (FWA) 
atn &A90A &A90A ЕМА = atn (FWA) 
cos &A98C &A990 FWA = cos (FWA) 
deg &ABEA &ABC5 FWA = deg (FWA) 
exp &AAB7 &AA94 FWA = exp (FWA) 
іп &A807 &A801 FWA = Ln CFWA) 
Log &ABDO &ABAB FWA = Log (FWA) 
pi &ABFO &ABCB FWA = PI 

rad &ABD9 &ABBA FWA = rad (FWA) 
sin &A997 8А998 FWA = sin (FWA) 
sqr &A7B7 &A7B7 FWA = sqr (FWA) 
tan &A6CC &A6C1 FWA tan (FWA) 
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5.2 Mathematical Functions Routines Description 


subroutine name : 
: FWA - acs (FWA) 

: &A8CF then &A929 
: &A8DD then &A927 


function 
BASIC 1 address 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


acs 


: FWA contains a floating point number between 


-1 and 41 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 32567 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: &A8DD 


BASIC 2 address 


entry conditions 


exit status 


typical timing 


: asn 


: FWA = asn (FWA) 


&A8CF 


: FWA contains a floating point number between 


—1 and +1 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 31970 microseconds 
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subroutine name 


function 


BASIC 2 address 


entry conditions 


exit status 


typical timing 


: atn 

: FWA = atn (FWA) 
BASIC 1 address : 
: &A90A 


&A90A 


: FWA contains a floating point number between 


-1Е38 and +1E38 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 19110 microseconds 
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subroutine name 


function 


BASIC 2 address 


entry conditions 


exit status 


error reports 


typical timing 


: cos 
: FWA = cos (FWA) 
BASIC 1 address : 
: &A990 


&A98C 


: FWA contains number of radians in floating 


point 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: accuracy lost 


: 26787 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: deg 

: FWA = deg (FWA) 
: &ABEA 

: «АВС5 


: FWA contains number of radians in floating 


point 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 1711 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: exp 
: FWA = exp (FWA) 

: &AAB7 

: «АА94 

: FWA contains a valid floating point argument 
: IWA unchanged 

: FWA = result 

: FWB destroyed 

: SWA unchanged 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: exp range 


: 14997 microseconds 
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subroutine name : 
: FWA = In (FWA) 
: &A807 
: &A801 


function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


In 


: FWA contains a valid floating point argument 
: IWA unchanged 

: FWA = result 

: FWB destroyed 

: SWA unchanged 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: log range 


: 17192 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: log 

: FWA = log (FWA) 
: &ABDO 

: &ABAB 

: FWA contains a valid floating point argument 
: IWA unchanged 

: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: log range 


: 1551 microseconds 
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subroutine name 


function 


BASIC 2 address 
entry conditions 
comments 


exit status 


typical timing 


: pi 

: FWA = PI 
BASIC 1 address : 
: &ABCB 


&ABFO 


: none 
: FWA set to 3.14159265 
: IWA unchanged 

: FWA = result 

: FWB unchanged 

: SWA unchanged 

: А destroyed 

: X unchanged 

: Y destroyed 

: P destroyed 

: &4B,&4C destroyed 


: 87 microseconds 
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subroutine name 


function 


BASIC 2 address 


entry conditions 


exit status 


typical timing 


: rad 

: FWA = rad (FWA) 
BASIC 1 address : 
: &ABB4 


&ABD9 


: FWA contains number of degrees in floating 


point 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 1566 microseconds 
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subroutine name 


function 


BASIC 2 address 


entry conditions 


exit status 


error reports 


typical timing 


: sin 

: FWA = sin (FWA) 
BASIC 1 address : 
: &A99B 


&A997 


: FWA contains number of radians in floating 


point 


: IWA unchanged 
: FWA = result 

: FWB destroyed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: accuracy lost 


: 15483 microseconds 
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subroutine name 


function 


BASIC 2 address 
entry conditions 


exit status 


error reports 


typical timing 


: sqr 
: FWA = sqr (FWA) 
BASIC 1 address : 
: &A7B7 


&A7B7 


: FWA contains a valid floating point argument 
: IWA unchanged 

: FWA = result 

: FWB destroyed 

: SWA unchanged 

: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 

: —ve root 


: 8783 microseconds 
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subroutine name : tan 

function : FWA = tan (FWA) 
BASIC 1 address : «А6СС 

BASIC 2 address : &A6C1 


entry conditions : FWA contains number of radians in floating 
point 


exit status : IWA unchanged 
: FWA = result 
: FWB destroyed 
: SWA unchanged 
: А destroyed 
: X destroyed 
: Y destroyed 
: P destroyed 
error reports : accuracy lost 


typical timing : 41770 microseconds 
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5.3 Mathematical Functions Demonstration 


The following program draws a circle in mode 0, centred at 600,500 
and with a radius of 400. It uses the polygon method, with 100 
sides in the figure. In BASIC, the program could be written as: 


10 MODE 0 
20 xcentre = 600 
50 ycentre = 500 


40 increment = 2xPI/100 

50 stop = 2xPI 

60 radius = 400 

70 MOVE xcentre+radius,ycentre 

80 FOR angle = 0 TO stop STEP increment 


90 DRAW xcentre + radius * COSCangle) , 
ycentre + radius ж SIN(angle) 

100 NEXT 

110 END 


It will be seen that the equivalent BASIC 2 assembly language code 
is somewhat more long-winded to write. Paradoxically, the 
machine code generated occupies less memory and runs a bit faster 
than the BASIC code above, in spite of the fact that it was written 
with clarity as the main objective rather than efficiency. 


Perhaps the most striking feature of the assembly language code is 
that it is all so simple, the clever code being all within the BASIC 
subroutines called. 


10 aunp = &A3B5:aplus = &A500 
20 apack = RA38D:atest = &9A5F 
30 amult = RA656:fpi1 = &A3E4 
40 cos = &A990:sin = &A99B 


50 oswrch = &FFEE 

60 DIM пс2 500 

70 FOR pass% = O TO 2 STEP 2 
80 P% = mc% 

90 COPT pass% 

100 .xcentre 


110 ОРТ FNEQUF(600) 
120 .xadd EQUW xcentre 

130 .ycentre 

140 ОРТ FNEQUF(500) 
150 .yadd EQUW ycentre 

160 .radius 

170 ОРТ FNEQUF(40D) 
180 .radadd EQUW radius 

190 .angle 

200 ОРТ  FNEQUF(0) 


210 .angadd EQUW angle 
220 .increment 
230 ОРТ FNEQUF(2*PI/100) 
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240 .incadd EQUW increment 
250 .stop 

260 ОРТ FNEQUF(2*PI) 
270 .stopadd EQUW stop 

280 .plot EQUB 25 

290 .рагт1 EQUB 4 

300 .xcoord  EQUW 1000 

310 .ycoord EQUW 500 

320 .тоде0 EQUW &16 


330 \ 

340 .start 

350 LDA modeO \ MODE 0 

360 JSR osurch \ 

370 LDA тоае0+1 \ 

380 JSR osurch \ 

390 JSR doplot \ MOVE 

400 LDA #5 \ reset for 

410 STA parmi \ PLOT 

420 .mainloop 

450 JSR pointincr N point &AB,&AC at increment 
440 JSR aunp \ unpack into FWA 

450 JSR pointangle N point &4B,&4C at angle 
460 JSR aplus \ add them 

470 JSR apack \ put result іп angle 
480 JSR pointstop N point &AB,&AC at stop 
490 JSR atest \ test for end 

500 BCC alldone \ yes - alldone 

510 JSR doxcoord \ calculate X coordinate 
520 JSR doycoord \ calculate Y coordinate 
530 JSR doplot \ plot it 

540 JMP mainloop \ and back 

550 .alldone 

560 RTS \ bye bye 

570 \ 

580 .doxcoord \ CALCULATE X COORD 

590 JSR pointangle \ point &4B,&4C at angle 
600 JSR aunp \ unpack into FWA 

610 JSR cos \ get cosine 

620 JSR pointrad N point &AB,&AC at radius 
630 JSR amult \ multiply 

640 JSR pointx \ point &4B,&A4C at xcentre 
650 JSR aplus \ add 

660 JSR fpi1 \ convert to integer 

670 LDX #0 \ set Loop counter 


680 .doxloop 


690 LDA &2A,X \ get next byte of IWA 
700 STA xcoord,X \ save it 

710 INX \ bump X 

720 CPX #4 \ are we done ? 

730 BCC dox Loop \ no – back 

740 RTS \ back 

750 \ 


123 


760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 


.doycoord 
JSR pointangle 
JSR aunp 
JSR sin 
JSR pointrad 
JSR amult 
JSR pointy 
JSR aplus 
JSR fpi1 
LDX #0 

-doyloop 
LDA &2A,X 
STA ycoord,X 
INX 
CPX #4 
BCC doy Loop 
RTS 

\ 

.doplot 
LDX #0 

.pLotLoop 
LDA plot,X 
JSR oswrch 
INX 
CPX #6 
BCC plotloop 
RTS 

\ 

.pointincr 
LDA incadd 
STA &4B 
LDA incadd+1 
STA &4C 
RTS 

\ 

-pointangle 
LDA angadd 
STA &4B 
LDA angadd+1 
STA &4C 
RTS 

\ 

.pointstop 
LDA stopadd 
STA &4B 
LDA stopadd+1 
STA &4C 
RTS 

\ 

-pointrad 
LDA radadd 
STA &4B 


TT Z Z Z Z Z Z Z 


C ыы ы” м ше 


7 


аа е ut utn ga Sa 


M rS QUERI Ық сей 


CALCULATE Y COORD 


point &4B,&4C 


at angle 


unpack into FWA 


get sine 
point &4B,&4C 
multiply 
point &4B,&4C 
add 


at radius 


at ycentre 


convert to integer 
set Loop counter 


get next byte 
save it 

bump X 

are we done ? 
no - back 
back 


PLOT 


of IWA 


zeroise loop conter 


get next byte 
print it 

bump X 

end of plot ? 
no - back 


POINT %4В,%4С 
Lo address of 
save 
hi address of 
save 
back 


POINT %4В,%4С 
Lo address of 
save 
hi address of 
save 


POINT %4В,%4С 
Lo address of 
save 
hi address of 
save 
back 


POINT &4B,84C 


Lo address of 
save 
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of plot 


at increment 
increment 


increment 


at angle 
angle 


angle 


at stop 
stop 


stop 


at radius 
radius 


1280 LDA radadd+1 
1290 STA &4C 
1300 RTS 

1310 \ 

1320 .pointx 

1330 LDA xadd 
1340 STA &4B 
1350 LDA xadd+1 
1360 STA &4C 
1370 RTS 

1380 \ 

1390 .pointy 

1400 LDA yadd 
1410 STA &4B 
1420 LDA yadd+1 
1430 STA &AC 
1440 RTS 

1450 1 

1460 NEXT pass% 

1470 CALL start 

1480 END 

1490 DEF FNEQUF(Z) 
1500 1% = 5 +?&&B4 +256*?&4B5 
1510 FOR JZ = 1 TO 5 
1520 ?P% = 21% 
1530 PZ = Р +1 
1540 IX = I4 +1 
1550 NEXT 


1560 = pass% 


hi address of 
save 
back 


POINT %4В,%4С 
lo address of 
save 
hi address of 
save 
back 


POINT %4В,%4С 
lo address of 
save 
hi address of 
save 
back 
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radius 


at xcentre 
xcentre 


xcentre 


at ycentre 
ycentre 


ycentre 


126 


6 RANDOM NUMBERS 


The BBC Micro generates random numbers by applying a pseudo- 
randomising algorithm to an initial value, known as the seed. After 
each application of the algorithm, the resulting random number is 
not only returned to the user, it also becomes the new seed. 
Naturally, facilities are also provided for the user to supply a value 
to initialise the seed. 


On power-up, the same seed is planted in the random number data 
field each time. For this reason, the algorithm cannot be completely 
random. The user can of course alter this by supplying a seed at the 
start of the program. Providing the user does not previously set 
TIME in the program, RND(-TIME) will achieve this. 


6.1 Random Number Work Area 


The area of memory dedicated to random numbers starts at &0D 
and ends at &11. This five byte area will be called the Random 
Number Work Area or RWA. It is the source data for all random 
number operations. In BASIC, there are many varieties of the RND 
statement. 


RND generates an integer random number between 
-2147483648 and --2147483647. It executes the 
algorithm and copies &00 to &10 into the IWA. 


RND(-X) resets the seed to X and returns -X. It copies the IWA 
into «ОП to &10, sets &11 to #&40 and leaves the IWA 
unchanged at -X. 


RND(0) repeats the last random number returned by RND(1). It 
simply copies the RWA into the FWA. 


RND(1) generates a floating point random number between 0 
and 0.999999. It executes the algorithm and copies the 
RWA into the FWA. 


RND(X) generates an integer random number between 1 and X. 
After executing the algorithm, the result is returned in 
the IWA. 
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6.2 Random Number Routine Summary 


The individual functions supported by the single BASIC statement 
RND have different entry points. Therefore, in assembly language 
programming, it is necessary to regard each function as a separate 
routine. À program is supplied which demonstrates each routine. 


Name BASIC 1 BASIC 2 Function 

address address 
гпа0 &AF9B &AF6C RND(0) repeat Last rnd1 
rnd1 &AF98 &AF69 RND(1) FWA = from О to 0.999999 
rndi &AF80 &AF51 RND IWA = random integer number 
rndseed &AF6E &AF3F RNDC-X) RWA seeded with X 
rndx &AF53 &AF24 RNDOO IWA = from 1 to X 
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6.3 Random Number Routines Description 


subroutine name 
function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


: гпао 

: FWA = value returned by last rnd1 
: &AF9B 

: &AF6C 

: none 

: IWA destroyed 

: FWA = result 

: FWB unchanged 
: RWA unchanged 
: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 82 microseconds 


129 


subroutine name : 


: FWA = random number from 0 to 0.999999 


function 

BASIC 1 address 
BASIC 2 address 
entry conditions 


exit status 


typical timing 


rnd1 


: &AF98 

: &AF69 

: попе 

: IWA destroyed 

: FWA = result 

: FWB unchanged 
: RWA changed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 793 microseconds 
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subroutine name : rndi 


: IWA = a random integer number from 
-2147483648 to +2147483647 


function 


BASIC 1 address : 
BASIC 2 address : 


entry conditions 


exit status 


typical timing 


&AF80 
&AF51 


: none 
: IWA = result 

: FWA unchanged 
: FWB unchanged 
: RWA changed 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 745 microseconds 
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subroutine name 
function 

BASIC 1 address 
BASIC 2 address 


entry conditions 


exit status 


typical timing 


: rndseed 

: RWA = IWA + #&40 in fifth byte 

: &AF6E 

: &AF3F 

: IWA set to value to be seeded. 

: Unlike in BASIC, it does not have to be 


negative. 


: IWA unchanged 
: FWA destroyed 
: FWB unchanged 
: RWA re-seeded 

: SWA unchanged 
: А destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 42 microseconds 
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subroutine name 


function 


BASIC 1 address : 
: &AF24 


BASIC 2 address 


entry conditions 


exit status 


typical timing 


: rndx 


: IWA = random number from I to value in IWA 


&AF53 


: IWA = maximum value of random number 
: &4,5 should be pointed to HIMEM 
: IWA - result 

: FWA destroyed 

: FWB unchanged 

: RWA changed 

: SWA unchanged 

: A destroyed 

: X destroyed 

: Y destroyed 

: P destroyed 


: 2974 microseconds 
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6.4 Random Number Demonstration 


The following program, written in BASIC 2, demonstrates all the 
random number functions. Note that the last routine to be tested is 
rndseed. Repeated runs of this program will produce identical 
results because of this. Delete lines 1030 to 1070 for more random 
results. Note also that гпао does indeed return the same result as 
the preceding rnd1. 


10 rndi = &AF51:rndseed = &AF3F 
20 rnd0 = &AF6C:rnd1 &AF69 
30 rndx - &AF24:apack - &A38D 
40 DIM mc% 200 

50 FOR pass% = 0 TO 2 STEP 2 

60 P% = mc% 

70 ГОРТ pass% 

80 .seedval  EQUD -123456 

90 .rndxval  EQUD 2345678 

100 .result EQUD O 


110 EQUB 0 

120 \ 

130 .rnditest 

140 JSR rndi \ do rndi 

150 JSR saveiwa \ save IWA in result 
170 RTS \ bye bye 

180 \ 

190 .rndseedtest 

200 LDX #0 \ zeroise Loop counter 


210 .rndseedloop 


220 LDA seedval,X \ get next byte of seed 
230 STA &2A,X \ save in IWA 

240 INX \ bump X 

250 CPX #4 \ end of Loop ? 

260 ВСС rndseedloop \ no - back 

270 JSR rndseed \ seed IWA into RWA 
280 RTS \ bye bye 

290 \ 

300 .rnditest 

310 JSR rnd1 \ do rnd1 

320 JSR savefwa \ save FWA in result 
330 RTS \ bye bye 

340 \ 

350 .rndOtest 

360 JSR rndO \ do rndO 

370 JSR savefwa \ save FWA in result 
380 RTS \ bye bye 

390 \ 

400 .rndxtest 

410 LDA 86 \ HIMEM Lo 

420 STA &4 \ save in &4 

450 LDA 87 \ HIMEM hi 

440 5ТА %5 \ save іп %5 
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450 LDX #0 
460 .гпахіоор 


ет 


zeroise loop counter 


470 LDA гпахуа(,Х \ get next byte of rndxval 
480 STA %2А,Х \ save in IWA 

490 INX \ bump X 

500 CPX #4 \ end of Loop ? 

510 ВСС гпах(оор \ no - back 

520 JSR rndx \ do rndx 

530 JSR saveiwa \ copy IWA into result 
540 RTS \ bye bye 

550 \ 

560 .saveiwa 

570 LDX 40 \ zeroise Loop counter 


580 .(ооріма 


590 LDA &2A,X \ get next byte of IWA 
600 STA result,X \ save in result 

610 INX \ bump X 

620 CPX #4 \ end of Loop ? 

630 ВСС Loopiwa \ no - back 

640 RTS \ back 

650 \ 

660 .savefwa 

670 LDA #result MOD 256 \ Lo address of result 
680 STA &4B \ save 1% 

690 LDA #result DIV 256 \ hi address of result 
700 STA &4C \ save 1% 

710 JSR apack \ pack FWA into result 
720 RTS \ back 

730 1 

740 NEXT passX 

750 CLS 

7602-20 

770 PRINT 


780 PRINT "Test of rnd1" 

790 CALL rnditest 

800 REM set 7 = result 

810 1% = 3 + ?&4B4 + 256*?&4B5 
820 FOR JZ = O TO 4 

830 ?(1%+J%) = ?(result+J4) 
840 NEXT 

850 PRINT "result returned = ";Z 
860 PRINT 

870 PRINT "Test of rnd0" 

880 CALL rndOtest 

890 REM set Z - result 

900 1% = 5 + ?&4B4 + 256*?&4B5 
910 FOR JZ = O TO 4 

920 2(IZ+JZ) = 2Cresult+JZ) 
930 NEXT 

940 PRINT "result returned = ";Z 
950 PRINT 

960 PRINT "Test of rndi" 
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970 CALL rnditest 


980 PRINT "result returned = 


990 PRINT 


1000 PRINT "Test of rndx" 


1010 CALL rndxtest 


1020 PRINT "result returned - 


1030 PRINT 


1040 PRINT "Test of rndseed" 
1050 CALL rndseedtest 

1060 PRINT "RWA bytes 0-5 = 
1070 PRINT "RWA byte 4 = 


1080 END 
>RUN 


Test of rnd1 
result returned 


Test of гпа0 
result returned 


Test of rndi 
result returned 


Test of rndx 
result returned 


Test of rndseed 
RWA bytes 0-3 
RWA byte 4 


0.424982422 


0.424982422 


"rlresult 


"rlresult 


"мер 
"2811 


-1.54400018Е9 


2241095 


-123456 
64 
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7 BASIC MEMORY МАР 


Descriptions of many of the memory areas used by the BASIC 
interpreter have already been given. This memory map is included 
to allow the reader to explore further into the BASIC ROM if 
desired. Doubtless there are other subroutines not mentioned in 
this book which the reader could use in specific applications. 


7.1 Zero Page Dedicated Locations 


These addresses are used identically in BASIC 1 and BASIC 2. АП 
two-byte address pointers are low,high. 


&00 - &01 LOMEM 
pointer to the start of BASIC variables. 


&02 - 805 VARTOP 
pointer to the end of BASIC variables. 


&04 - &05 BASIC STACK POINTER 


pointer to most recent entry in the BASIC stack. 


&0A BASIC TEXT POINTER OFFSET 


the offset with respect to &B,&C of the 
byte of BASIC text currently being processed. 


&0B - &0C BASIC TEXT POINTER 


pointer to start of BASIC text Line being processed. 


&0D — 811 RND WORK AREA 
see Chapter 6. 


pointer to the end of BASIC program not 
including variables. 
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&14 PRINT BYTES 
the number of bytes іп a print output field. 


&15 PRINT FLAG 


0 = decimal 
-ve = hexadecimal 


816 - 817 ERROR ROUTINE VECTOR 


pointer to the address of the BASIC error routine. 


&18 PAGE DIV 256 


page number where BASIC program starts. 


&19 - &1A SECONDARY BASIC TEXT POINTER 
secondary &В,&С. 


&1B SECONDARY BASIC TEXT OFFSET 


secondary &A. 


&1C - &1D BASIC PROGRAM START 
pointer to start of BASIC program. 


&1F LISTO FLAG 
number ANDed from the List below: 


a 
0 - LISTO off 

1 = insert space after Line number 
2 - indent FOR Loops 

4 = indent REPEAT Loops 


&20 TRACE FLAG 


0 = trace off 
1 = trace on 


&24 REPEAT LEVEL 
number of nested REPEATs outstanding. 
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825 GOSUB LEVEL 


number of nested GOSUBs outstanding. 


&26 15*FOR LEVEL 


15 * number of nested FOR Loops outstanding. 


&27 VARIABLE TYPE 


&00 - byte 

&04 - integer 

&05 - floating point 
&81 - string 

&A4 - function name 
&F2 - procedure name 


bit 0- List flag 
bit 1 = errors flag 
bit 2 = relocate flag (BASIC 2 only) 


7.2 Zero Page Multiple Use Locations 


These addresses are used identically in BASIC 1 and BASIC 2. The 
main uses only are given. 
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7.3 Resident Integer Variables 


АП are stored with the least significant byte first. 


&400 - &403 ах 
&404 - &407 АЖ 
&408 - 840B B% 
&40C - &40F CX 
&410 - &413 DZ% 
&414 – &417 EX 
8418 - &41B F% 
&41C — &41F GX 
&420 - &423 H% 
9424 — 8427 1% 
9428 - &42B J% 
&42C - &&2F КЖ 
8430 – 9455 L% 
9454 - 8457 MX 
98458 - &43B NX 
&43C - 843F 0% 
8440 - 9445 P% 
9444 - 8447 02 
8448 — &&&B ЮЖ 
&44C — &&&F S% 
8450 - &453 T% 
&454 — 8457 U% 
&458 - &45B VX 
&45C — &45F W% 
&460 - &465 XX 
&464 - 8467 Y% 
8468 - &46B 72 


7.4 Floating Point Temporary Areas 


АП are stored in 5 byte packed floating point format. 


&46C - &470 TEMP 1 
&471 — %475 TEMP 2 
&476 - &47A TEMP 3 
&47B - &47F TEMP 4 
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7.5 Variable Pointer Table 


There is a variable look-up table for each character with which a 
variable name can start. Each pair of addresses is a lo,hi pointer to 
the variables whose names start with a particular character. 


9480-8481 = à &4AA-&4AB = U &4D2-84D3 = i 
9482-8485 = A &4AC-&4AD = V &ADA-&AD5 = j 
9484-8485 = B &4AE-&4AF = W &4D6-&4D7 = k 
9486-8487 = С &4B0-&4B1 = X 9408-8409 = L 
9488-8489 = D &AB2-&ABS = Y &ADA-&ADB = m 
&48A-848B = E &4B4-84B5 = Z &4DC-&4DD = n 
&48C-&48D = F &4B6-&4B7 = Г &ADE-&ADF = о 
&48E-&48F = G &4B8-&4B9 = \ &4E0-&4E1 = p 
9490-8491 = H &ABA-&ABB = 1 &AE2-&AES = q 
9492-8495 = І &4BC-&4BD = ^ &4E4-84E5 = r 
9494-8495 = J &4BE-&4BF = _ &4E6-&4E7 = s 
9496-8497 = K &4CO-R4C1 = £ 94Е8-84Е9 = t 
9498-8499 = L &4C2-R4C3 = а &4EA-&4EB = u 
&49A-849B = М &4C4-84C5 = b &4EC-&4ED = v 
&49C-&A9D = М &4C6-&4C7 = c &4EE-&4EF = w 
&A9E-&A9F = 0 &4C8-&AC9 = d &AFO0-&AF1 = x 
&4AO-R4A1 = P &ACA-&ACB = e &4F2-84F3 = y 
94А2-%4АЗ = 0 &ACC-&ACD = f &AFA-RAF5 = z 
&4A4-84A5 = R &ACE-&ACF = g &4F6-&4F7 = procedures 
&4A6-&4A7 = S &4DO-&4D1 = h &4F8-&4F9 = functions 
94А8-%4А9 = Т 


The use of these look-up tables сап best be explained by ап 
example. Suppose we run the following program: 


ТО DIM АККАҮ2(3) 

20 FOR AZ = 0 TO 3 
30 ARRAYZ(AZ) = AX 
40 NEXT 

50 ALPHA$ = "TESTING" 
60 AFPNUM = 1.0 

70 AINTZ = 1 

80 END 


Assume that &482 contains #&0E and &483 contains #&16. Then 
the look-up table for variables starting with the letter ‘A’ will be 
found at &160E which will contain the following: 


&160E - Я%2А pointer to next variable 
&160F = #&16 at &162A 

&1610 = 7852 R rest 

&1611 = #952 R of 

&1612 = #&41 А name 

&1613 = #859 Y for 

81614 = H&25 X  ARRAYZ( 

&1615 = #&28 ( 
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81616 = Я%00 end of name marker. 


81617 = #803 2 * number of dimensions + 1 
81618 = Я%04 number of elements іп 15% dimension 
81619 = Я%00 number of elements іп 2nd dimension 
&161A = Я%00 contents 

&161B = #&00 of 

&161C = #&00 ARRAYZ(0) 

&161D = #&00 = 0 

&161E = #801 contents 

&161F = #&00 of 

&1620 = #&00 ARRAYZ%(1) 

&1621 = #&00 = 1 

81622 = #&02 contents 

81605 = #&00 of 

81604 = #800 ARRAY%(2) 

&1625 = #&00 =2 

81626 = Я%05 contents 

&1627 = #&00 of 

81608 = H&00 ARRAY%(3) 

&1629 = #800 = 3 

&162A = #&3D pointer to next variable 

&162B = #&16 at &163D 

&162C = АС L rest 

9160) = #850 Р of 

&162E = #848 Н variable 

&162F = #&41 А name 

&1630 = #824 $ for ALPHA$ 

81651 = Я%00 end of name marker 

81652 = #836 pointer to 81656 which contains 
81655 = H&16 current value of ALPHA$ 

81654 = Я%07 maximum size allocated to ALPHA$ 
&1635 = Я%07 current size of ALPHA$ 

81656 = #&54 Т contents 

91637 = #845 E of 

&1638 = #853 5  ALPHA$ 

91639 = Я854 T 

&163A = H&49 I 

&163B = #&4E N 

&163C = #&47 G 

8165) = #&4A pointer to next variable 

&163E = #&16 at &164A 
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&163F = #846 F rest 

&1640 = #&50 Р of 

&1641 = #&4E N name 

&1642 = #&55 U for 

&1643 = #&4D M  AFPNUM 

81644 = Я%00 end of name marker 
81645 = #&81 contents of 

81646 = Я%00 AFPNUM 

91647 = #&00 = +1.0 

81648 = Я%00 іп 5 byte, packed, 
81649 - Я%00 floating point format 
&164A = H&OO pointer to next variable - 80000 
&164B - Я%00 so Last entry іп table 
&164C = #849 I rest of 

&164D = #&84E М variable 

&164E = #854 Т name for 

8164F = #825 2 АТАТ? 

81650 = Я%00 end of name marker 
&1651 = #&01 AINTZ 

&1652 - Я%00 contains 

&1653 = #&00 +1 in 4byte, 

&1654 - Я%00 integer format 


Each entry in the look-up table starts with a 2 byte (lo,hi) pointer to 
the next entry in the table, except for the last entry which contains 
#&0000. Following this is the rest of the variable name, excluding 
its initial letter, and then a zero byte to indicate the end of the 
name. This is followed by a contents section, the exact format of 
which depends on the type of variable being stored. Thus integer 
variables are stored in 4 bytes, whilst floating point variables are 
stored in 5 bytes. String variables can change in size during the 
running of a program. Thus the table contains a pointer to the 
string (rather than the string itself), together with details of its 
allocated and current size. Arrays can be integer, floating point or 
string. Integer and floating point array contents are stored in the 
table, whilst string arrays point to the string. At the start of the 
contents, there are some additional bytes which define the size of 
the array. 
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7.6 BASIC Stacks апа Buffers 


&500 - &5FF FOR/REPEAT/GOSUB STACK 
&600 - &6FF STRING WORK AREA / CALL PARAMETER BLOCK 
8700 - &7FF BASIC LINE INPUT BUFFER 


BASIC reads text entered at the screen into the line input buffer. 
The text is then tokenised (each BASIC command is replaced by a 
number). If the input line starts with a line number, the tokenised 
line is inserted into a BASIC program at the proper place. 
Otherwise the tokenised line is executed immediately. 


7.7 BASIC Token and Action Tables 


Tokens range from &80 to &FF. BASIC has a token table which for 
each BASIC command contains: 


a) the command name in ASCII 
b) the token 
c) а flag byte used by the interpreter 


The token table is located at: 


&806D - 88558 BASIC 1 
88071 - 8856С BASIC 2 


The order of commands in the table is important since it specifies 
the minimum acceptable abbreviation of a command. À command 
will be recognised so long as it starts with one or more letters of 
that command and ends in a full-stop, provided the abbreviation 
does not match a command earlier in the table. The first command 
in the table is ‘AND’ for which ‘A.’ is sufficient. The second is 
‘ABS’. ‘A.’ cannot be used for ABS because the previous command, 
‘AND’, is matched by it. However ‘AB.’ is quite satisfactory. 


&8D is a special token used to prefix a BASIC line number. 
Commands lower than &8E (&8F in BASIC 1) are processed by the 
interpreter as and when they occur in a line. The rest of the 
commands have an associated action address found by looking up a 
two-part action address table. The first part contains the low byte 
addresses, whilst the second contains the high bytes. Each table is 
indexed by the token number less &8E (&8F in BASIC 1). 


The action address tables are located at: 


Lo bytes &835A - 885СА BASIC 1 
8856) - &83DE BASIC 2 


hi bytes &83CB - 8845С BASIC 1 
&83DF - 88450 BASIC 2 


The following tables summarise BASIC commands іп alphabetic 
order. 
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7.8 BASIC Tables Summary 


BASIC 
COMMAND 


CLS 
COLOUR 
cos 
=COUNT 
DATA 
DEF 
DEG 
DELETE 
DIM 
DIV 
DRAW 


BASIC 1 BASIC 2 TOKEN 
ADDRESS ADDRESS НЕХ 
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ВА5ІС 
COMMAND 


ENDPROC 
ENVELOPE 
EOF 
EOR 
ERL 
ERR 
ERROR 
EVAL 
EXP 
EXT 
FALSE 
FN 

FOR 
GCOL 
GET 
GETS 
GOSUB 
GOTO 
=HIMEM 
HIMEM= 
IF 
INKEY 
ІМКЕҮ% 
INPUT 
INSTRC 
INT 


BASIC 1 BASIC 2 TOKEN 
ADDRESS ADDRESS НЕХ 
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ВА5ІС 
COMMAND 


LEFT$( 
LEN 
LET 
LINE 


Line no. 


LIST 
LN 
LOAD 
LOCAL 
106 
=LOMEM 
LOMEM- 
МІр%( 
MOD 
MODE 
MOVE 
NEW 
NEXT 
NOT 
OFF 
OLD 

ON 
OPENIN 
OPENOUT 
OPENUP 
OR 


BASIC 1 BASIC 2 TOKEN 
ADDRESS ADDRESS НЕХ 


----- &BF80 &AD 
----- --- 884 
----- &BEC2  &FF 


&AD in BASIC 1 
BASIC 2 only 


BASIC 2 only 
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ВА5ІС BASIC 1 BASIC 2 ТОКЕМ 
COMMAND ADDRESS ADDRESS НЕХ 


PAGE= 89239 89283 &DO 
PI &ABFO  &ABCB &AF 
PLOT &93AE &93F1 &FO 
РОІМТ( &AB64  RAB41 %В0 
POS &AB92  &AB6D &B1 
PRINT 88035 %809А &F1 
PROC &92B6 89304 &F2 
=PTR &BF50 &BF47 &8F 
PTR= &BF39 &BF30 &CF 
RAD &ABD6  RABB1 &B2 
READ &BB39 &BB1F &F3 
REM &8AED &8B7D &F4 


RENUMBER &8F37 &8FA3 &CC 
REPEAT &BBFF  &BBE4 &F5 
REPORT &BFEÓ &BFE4 &F6 
RESTORE %8В00  &BAE6 &F7 
RETURN &B8D5 &B8B6 &F8 
RIGHT$( &BO1D &АҒЕЕ &C2 
RND &AF78  &AF49 &B3 
RUN &BD29  &BD11 &F9 
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ВА5ІС 
COMMAND 


STRINGS 
ТАВ( 
ТАМ 
ТНЕМ 
-ТІМЕ 
ТІМЕ- 
TO 
TRACE 
TRUE 
UNTIL 
USR 


BASIC 1 BASIC 2 TOKEN 
ADDRESS ADDRESS НЕХ 
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8 TIMINGS 


Frequently the speed of a program determines its success or failure, 
especially in graphics applications. Thus, the execution time of a 
program is an important topic, not least to those who decide to use 
the technique advocated in this book, and a separate chapter on 
this topic is justified. 


8.1 Units of time 


The following units of time are frequently encountered when 
dealing with computers: 


1 picosecond = ----------------- seconds 
1,000,000,000,000 
1 
1 nanosecond = ------------- seconds 
1,000,000,000 
1 
1 microsecond = --------- seconds 
1,000,000 
1 
1 millisecond = —-—— seconds 
1,000 
1 
1 centisecond = --- seconds 
100 


8.2 Computer Processor Speed 


The speed of a cpu is determined by its cycle time. A cycle is the 
interval of time that elapses between successive pulses of the 
system clock. During a cycle the cpu performs a fundamental cpu 
operation. This operation could be to fetch a byte of data from RAM 
or to store a byte of data in RAM, for example. Each machine code 
instruction takes several cycles. The more complicated instructions 
take more cycles than the simpler ones. In fact the number of cycles 
needed for a 6502 instruction is always between 2 and 8. 


The BBC micro has a cycle time of 0.5 microseconds, which is just 
another way of saying that it is clocked at 2 Megahertz (2 million 
times per second). Thus it can perform its fastest instructions 

(2 cycles) іп 1 microsecond, whilst its slowest instructions 

(8 cycles) take 4 microseconds. 
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8.3 Program Speed 


Because of the importance of program speed, typical timings have 
been given for each of BASIC's subroutines. Note that the time 
spent in a given subroutine is not a constant. The exact number of 
instructions executed will vary a little depending on the data 
values being acted upon. Nevertheless, these timings can be used to 
forecast whether or not a program is going to be fast enough 
without actually writing the program. 


Thus, in the polygonal circle in chapter 5 it was recognised that in 
addition to some other code, the program would involve 100 sines 
and cosines. Timings given for sine and cosine predict that these 
functions alone will take 4.2 seconds. Broad-brush estimating such 
as this may well be sufficient to discard a technique. The 
calculation can be refined to give a more precise estimate if 
required. 


Moreover, timings given in this book can often be used to estimate 
run times of purely BASIC programs. In general, calling BASIC's 
subroutines from machine code saves about 1096 of the execution 
time. Thus by adding 1096 to the figures in this book it is possible 
to get a good approximation of the run time of many BASIC 
program commands also. 


8.4 Microsecond Timer 


Timings given in this book exclude any set up code required by the 
subroutine, but include time taken by interrupt routines such as 
servicing the various clocks maintained by the software. The BASIC 
TIME facility has a resolution of 1 centisecond which is 
insufficiently precise for many purposes. The following code uses 
the User 6522 to obtain 1 microsecond resolution of time and may 
be used to time another piece of code. The code to be timed is 
placed at ‘test’. Any set up code required, but not to be timed is 
placed at ‘setup’. 
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10 DIM mcZ 500 


20 FOR pass% = 0 TO 2 STEP 2 


30 PZ = тс# 

40 СОРТ pass% 

50 .save EQUW 0 
60 .time EQUD О 
70 .timer 

80 LDA 8206 

90 STA save 
100 LDA 8207 
110 STA save+1 
120 JSR setup 
130 SEI 


140 LDA #interrupt MOD 256 


150 STA &206 


160 LDA #interrupt DIV 256 


170 STA &207 
180 CLI 

190 LDA #&DF 
200 AND &FE6B 
210 STA &FE6B 
220 LDA #&A0 
230 STA &FEÓD 
240 STA &FEÓE 
250  JSR newinterrupt 
260  JSR test 
270 LDA &FEÓ8 
280 STA time 
290  LDA &FE69 
300 STA time+1 
310 SEI 

320 LDA save 
330 STA &206 
340 LDA save+1 
350 STA &207 
360 LDA #&20 
370 STA &FEÓD 
380 STA &FEÓE 


390  CLI 

400  RTS 

410 .interrupt 
420  CLC 

430 LDA time+2 
440 ADC #1 


450 STA time+2 

460 ВСС interruptdone 
470 LDA time+3 

480  ADC #0 

490 5ТА time+3 

500 .interruptdone 

510 JSR newinterrupt 
520  JMP (save) 


Y uto ut Gui Sao Gm aae ыса” э M wa cut ut 


P m т, AO uat pe ын” 


ALA Z Z 


save interrupt 

handler address Lo 

and save 

hi address too. 

set up. 

prevent interrupts. 
re-direct IRQ2V to 

our own 

routine called 
"interrupt". 

allow interrupts again. 
set up ACR of user 6522 
to put TIMER 2 

in count down mode. 
enable TIMER 2 
interrupts. 


reset TIMER 2. 

do code to be timed. 
get count down value 
and save in "time". 
also get 

hi byte 

prevent interrupts. 
restore IRQ2V. 


clear 

TIMER 2 interrupts. 
enable interrupts 

bye bye 

bump time by 1 

every time count 

down 

expires 

reset count down timer 


back to IRG2V. 
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530 .newinterrupt 


540 LDA #RFD \ reset 

550 STA &FEÓ8 \ TIMER 2 
560 LDA #&ЕЕ \ as #&FFFD 
570 STA &FE69 

580 RTS \ out 

1000 .test 

1010 \ CODE TO BE TIMED 

1020 RTS 

2000 .setup 


2010 \ SET UP FOR CODE TO BE TIMED 
2020 RTS 

3000 1:МЕХТ pass% 

3010 CALL timer 

5020 ?Ctime) = &FF - ?(time) 

3030 ?(time+1) = &FF - 2(time+1) 
3040 PRINT !time;" microseconds" 
3050 END 
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8.5 BASIC Timings 


The following tables list all 69 BASIC routines in alphabetic order, 
together with their typical timings. 


aclear š 25 
acomp Н 54 
acopyb 2 36 
acs : 32,567 
adiv - 1,545 
adiv10 : 360 
aminus š 254 
amult : 1,581 
amult1 š 1,508 
amult10 Н 171 
апогт : 27 
aone : 37 
араск 2 46 
apack1 Ч 51 
apack2 2 53 
apack3 : 55 
aplus : 246 
aplusb : 58 
aplus1 : 41 
arecip š 1,619 
around : 22 
ascnum : 1,748 
asign š 25 
asn š 31,970 
atest : 70 
atn : 19,110 
aunp : 49 
aunp1 : 62 
bclear : 25 
bcopya : 56 
bunp : 51 
cos š 26,787 
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deg 2 1,711 
ехр : 14,997 
fpascdec : 4,878 
fpaschex : 582 
fpi1 š 405 
fpi2 : 587 
iascdec š 5,369 
1а5сһех : 216 
icomp š 31 
idiv : 794 
ifpa : 149 
iin : 54 
iminus : 52 
imod Н 782 
imult Н 164 
ineg1 : 20 
jout : 57 
iplus : 50 
ipos : 17 / 27 
ismall : 20 
itest : 58 
izero š 25 
ігріп 27 
izpout 26 
іп 17,192 
[од 1,551 
pi 87 
rad 1,566 
rndO 82 
гпа1 2 795 
rndi : 745 
rndseed : 42 
гпах : 2,974 
sin Н 15,483 
sqr 8,783 
tan 41,770 
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9 TRIGONOMETRICAL 
MANIPULATIONS 


The previous chapter gave typical timings for all of the BASIC 
subroutines. Inspection of these timings reveals the fact that 
trigonometrical functions are especially time consuming. There are 
a number of different methods which can often be used to get round 
this problem and this chapter explains many of them. Each method 
is illustrated by the polygonal circle discussed in Chapter 5, but the 
methods are applicable to many situations in which trigonometry is 
used. It should be remembered that the conventional method of 
drawing circles takes nearly 6 seconds in BASIC and even in 
machine code requires 5.5 seconds. With a little chicanery, 
considerable improvements on these times may be achieved. Apart 
from the first method which must use assembly language, BASIC is 
used in demonstration programs so that the methods are easier to 
understand. 


9.1 Fixed Shapes Method 


The following program illustrates a method that can be used 
whenever the application draws a geometric shape of fixed 
dimensions in a fixed position. In this method, all the coordinates 
to be plotted are stored as constants within the program. In the 
demonstration program, the two functions, XCOORD and YCOORD 
respectively, store away all 100 Х and Ү coordinates of the circle to 
be plotted. The program itself simply plots these points. АП of the 
BASIC parts of the program are disposable. The generated machine 
code draws the circle in 0.28 seconds. Of course, this method uses a 
lot of memory to store coordinates (404 bytes in this example). 
Moreover, it is not a general purpose routine to draw many circles 
of different sizes. Nevertheless, it is the quickest method and is 
useful in many applications. 
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Fixed Circle Example 


0 MODE O:oswrch = &FFEE 

10 DIM mc% 1000 

20 FOR pass% = 0 TO 2 STEP 2 
30 PX = mc% 

40 СОРТ pass% 

50 .xcoords 


60 OPT FNXCOORD 
70 .ycoords 
80 OPT FNYCOORD 


90 .(оорсп% EQUB -1 
100 .plot EQUB 25 
110 .parms  EQUB 4 


120 EQUD O 

130 .circle 

140 LDA Z&16 \ set mode 

150 JSR osurch \ = 

160 LDA #0 \ zero 

170 JSR osurch 

180 .(оор 

190 LDX Loopcnt \ get Loopcnt 
200 INX \ bump it 

210 СРХ #101 \ test for 101 
220 BCS out \ if > 100 out 
230 5ТХ Loopcnt \ save Loopent 
240 TXA \ put X in A 
250 ASL A \ х 2 

260 ТАХ \ back іп Х 

270 LDA xcoords,X АХ get xcoord (о 
280 STA parms+1 

290 LDA ycoords,X N get ycoord (о 
300 STA parms+3 

310 INX \ bump X 

320 LDA xcoords,X АХ get xcoord hi 
330 STA parms+2 

340 ША ycoords,X N get ycoord hi 
350 STA parms+4 

360 LDX #0 \ reset plot Loop 
570 .plotloop 

380 LDA plot,X \ next byte of PLOT 
390 JSR osurch \ PLOT it 

400 INX \ bump X 

410 CPX #6 \ end of plot 
420 BCC plotloop \ no - back 
430 LDA #5 \ reset for PLOT 
440 STA parms 

450 JMP [оор \ next PLOT 

460 .out 

470 RTS \ bye bye 

480 J 


490 NEXT pass% 
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500 ТІМЕ-0 

510 CALL circle 

520 T = TIME / 100 

530 PRINT 

540 PRINT "time taken = ";T;" seconds" 
550 END 

560 DEF FNXCOORD 

570 FOR JZ- 0 TO 100 

580 XZ=INT (600+400*C0S(J%*2*PI/100)) 
590 ?P% = 28460 


600 P% = PX + 1 
610 2P% = 28461 
620 P% = РХ +1 
630 NEXT 


640 = pass% 

650 DEF FNYCOORD 

660 FOR JZ = 0 TO 100 

670 YZ=INT(500+400*SINCJ%Z*2*P1/100)) 
680 ?PZ = 28464 


690 PZ = РХ +1 
700 ЭРХ = 28465 
710 P% = РХ +1 
720 NEXT 


730 = pass% 


9.2 Reduced Accuracy Method 


In this method the number of sides in the polygon is reduced. The 
following BASIC program has only 32 sides in the polygon. The 
circle so drawn is quite adequate and is completed in 1.99 seconds: 


10 TIME = 0 

20 MODE 0 

50 хсепіге - 600 

40 ycentre = 500 

50 increment = 2xP1/32 

60 stop = 2xPI 

70 radius = 400 

80 MOVE xcentre+radius,ycentre 

90 FOR angle = 0 TO stop STEP increment 


100 DRAW xcentre + radius * COSCangle) , 
ycentre + radius ж SIN(angle) 

110 NEXT 

120 T=TIME/100 

130 PRINT 

140 PRINT "time taken = ";T;" seconds" 

150 END 
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9.3 Mathematical Transform Method 


In this method, the mathematics of the application are transformed 
in order to remove the most time-consuming mathematical 
functions from big loops. The following trigonometrical identities 
frequently prove useful: 


SINCA+B) = SINCA)*COS(B) + COSCA)*SIN(B) 
SINCA-B) = SINCA)*COS(B) - COSCA)*SIN(B) 
COSCA+B) = COSCA)*COS(B) - SINCA)*SIN(B) 
COSCA-B) = COSCA)*COS(B) + SINCA)*SIN(B) 


In the polygonal circle, the radius rotates like the sweep of a radar 
screen, the angle being incremented by 2*PI/N radians each time 
(N = number of sides in the polygon). Let us assume that at any 
given instant the sweeping radius is at an angle A radians. Let us 
also assume that ‘sold’ is the sine of A and ‘cold’ is the cosine of A 
at this moment. After one more increment, the radius will have 
swept through A+2*PI/N radians, which we will assume has a sine 
of ‘snew’ and a cosine of ‘cnew’. But from the identities above: 


SINCA+2*PI/N) = SINCA)*COS(2*PI/N) + COSCA)*SINC2*PI/N) 


or 

snew = soldxCOS(2xPI/N) + cold*SIN(2*PI/N) 
similarly 

cnew = coldxCOS(2xPI/N) - sold*SIN(2*PI/N) 


It is only necessary to calculate the sine and cosine of the 
increment, 2*PI/N, which can be done outside the big loop, since 
cnew and snew can then both be derived without recourse to 
further trigonometry. 


The following program, which draws a 100-sided circle in 2.48 
seconds, uses this technique. Note that at the start of the circle 
when the angle swept out is zero radians, the sine of zero is zero, 
but the cosine of zero is one (lines 90 and 100). The program is 
written in BASIC so that the transformation is clear. In assembly 
language about a further quarter of a second can be saved without 
any special effort. 


10 MODE 0 

20 ТІМЕ - 0 

30 xcentreX = 600 
40 ycentre% = 500 


50 radius% = 400 
60 sides% = 100 


70 sinc = SIN(2*PI/sides%) 
80 cinc = COS(2*PI/sides%) 
90 sold = 0 
100 cold = 1 
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110 MOVE xcentreZ*radiusAZ,ycentreZ 
120 FOR 17 = 0 TO sides% 


130 snew = soldxcinc + coldxsinc 
140 cnew = coldxcinc - soldxsinc 
150 DRAW xcentre% + radius% * cnew , 


ycentre% + radius% * snew 
160 sold = snew 
170 cold = cnew 
180 NEXT 
190 T=TIME/100 
200 PRINT : PRINT "time taken = ";T;" seconds" 
210 END 


9.4 Symmetry Method 


This method takes advantage of the symmetry of many geometric 
shapes. More than one point can be plotted on the basis of a single 
calculation. The following identities are used: 


SINCPI/2 + A) = + COS(A) 
SINCPI + A) -- SINCA) 
SINCS«PI/2 + A) = - COS(A) 
SINC2xPI + A) = + SIN(A) 
SINCPI/2 - A) = + COS(A) 
SINCPI - А) = + SINCA) 
SINC3xPI/2 - A) = - COS(A) 
SINC2xPI - A) = - SIN(A) 
COSCPI/2 + A) = - SIN(A) 
COS(PI + A) = – COSCA) 
COSC3xPI/2 + A) = + SIN(A) 
COSC2xPI + A) = + COS(A) 
COSCPI/2 - A) = + SIN(A) 
COS(PI — A) -- COSCA) 
COSC3xPI/2 - А) = - SIN(A) 
COSC2xPI - A) = + COS(A) 
TANCPI/2 + A) = - 1/TANCA) 
ТАМ(РІ + А) = + TAN(A) 
TAN(3*PI/2 + A) = - 1/TANCA) 
TAN(2*PI + A) = + ТАМА) 
TANCPI/2 - A) = + 1/TANCA) 
TANCPI - А) = - TANCA) 
TAN(3*PI/2 - А) = + 1/TANCA) 
ТАМС2*РТ - A) = - ТАМА) 


Thus the polygonal circle сап be plotted four points at a time, for 
example. Ав the radius sweeps round, not only the point at the 
current angle, but also points further round by 90 degrees (РІ/2 
radians), 180 degrees (PI radians) and 270 degrees (3*PI/2 radians) 
can be plotted without further trigonometry. 
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At any given time, if the angle swept out by the radius is A radians, 
the four coordinates will be: 


X1 = xcentreZ + radius#xCOS(A) 
Y1 = ycentre% + radius4xSIN(A) 


X2 = хсепігей + radius%*COS(PI/2+A) 
Y2 = усепігей + radiusZ*SIN(PI/2+A) 


ХЗ = xcentreX + radius%*COS(PI+A) 
ҮЗ = усепігей + radiusZ*SIN(PI+A) 


X4 = хсепігей + radius#*xCOS(3xP1/2+A) 
Y4 = усепігей + radiusZ*SIN(3*PI/2+A) 


The identities above give the following: 


X1 = xcentre% + radius%*COS(A) 
Y1 = ycentre% + radiusZ*SIN(A) 


X2 = xcentreX - radius4xSIN(A) 
Y2 = ycentre% + radius%*COS(A) 


ХЗ = хсепігей - radius%*COS(A) 
ҮЗ = ycentreX - radius4xSIN(A) 


X4 = xcentre% + radius#xSIN(A) 
Y4 = усепігей - radius%*COS(A) 


Thus it can be seen that four points can be plotted for a single value 
of A, knowing merely the sine and cosine of A. The four X 
coordinates are stored in an array called X, while the four Y 
coordinates are stored in an array called Y. Lines 70 to 160 in this 
program initialise both arrays with values corresponding to an 
angle of zero radians. The BASIC program, plotting four points at a 
time, draws a 100 sided circle in 3.11 seconds: 


10 MODE 0 

20 ТІМЕ - 0 

30 xcentre% = 600 

40 ycentre% = 500 

50 radius% = 400 

60 sides% = 100 

70 DIM X(4) 

80 DIM YC4) 

90 X(1) = xcentre% + radius% 

100 X(2) = xcentre% 
110 X(3) = xcentre% - radius% 
120 X(4) = xcentre% 
130 Y(1) = ycentreZ 
140 Y(2) = ycentre% + radius% 
150 YC3) = ycentre% 
160 Y(4) = ycentre% - radius% 
170 FOR 17 = 0 TO sides% DIV 4 
180 sinA = SIN(2*PI*I%Z/sidesz) 
190 cosA = COS(2*PI*I%/sidesz) 
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200 MOVE X(1),Y(1) 

210 X(1) = xcentre%+radius%*cosA 
220 Y(1) = ycentreZ*radiusX*sinA 
230 DRAW X(1),Y(1) 

240 MOVE X(2),Y(2) 

250 Х(2) = xcentreZ-radiusAX*sinA 
260 Ү(2) = ycentre%+radius%*cosA 
270 DRAW X(2),Y(2) 

280 MOVE X(3),Y(3) 

290 X(3) = xcentreZ-radius%*cosA 
300 Y(3) = ycentreX-radius#*xsinA 
310 DRAW X(3),Y(3) 

320 МОМЕ X(4),Y(4) 

330 X(4) = xcentreZ*radiusX*sinA 
340 Y(4) = ycentreZ-radiusX*cosA 
350 DRAW X(4),Y(4) 

360 NEXT 

370 T = TIME/100 

380 PRINT 

390 PRINT "time taken = ";T;" seconds" 
400 END 


This can be further improved. The circle can be drawn in 2.08 
seconds by plotting eight points at a time for each of the following 
angles: 


A radians 

PI/2 - A radians 
PI/2 * A radians 
PI A radians 
PI +A radians 
3xP1/2 - А radians 
3*PI/2 + A radians 
2xPI -А radians 


Note that the number of sides is increased to 104 simply to make it 
divisible by 8. Note also that in this example, the eight X and Y 
coordinates are not saved in an array, but are re-worked each time. 
This is simply an alternative method. 


The eight sets of coordinates to be plotted are: 


angle A X1 = хсепігей + radius%*COS(A) 
Y1 = ycentre% + radiusZ*SIN(A) 


X2 = хсепігей + radius#xSIN(A) 
Y2 = ycentre% + radius%*COS(A) 


angle PI/2 


І 
> 


angle РІ/2 + A ХЗ = хсепігей - radius4xSIN(A) 
ҮЗ = ycentre% + radius%*COS(A) 


angle PI -А X4 = хсепігей - radius%*COS(A) 
Y4 = ycentre% + radiusZ*SIN(A) 


angle PI + А X5 = хсепігей - radius%*COS(A) 
Y5 = усепігей - radius4xSIN(A) 
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І 
> 
> 
о 

I 


= xcentreX 
Y6 = ycentrex 


angle 3xP1/2 


angle 3*PI/2 + A X7 = хсепігей 
Y? = ycentrex 


angle 2*PI - А X8 = xcentrex 
Y8 = ycentrex 


- radiusX*SINCA) 
- radiusX*COS(CA) 


+ radiusX*SINCA) 
- radiusX*COS(CA) 


+ radius%Z*COS(A) 
- radiusX*SINCA) 


10 ТІМЕ-0 

20 MODE 0 

30 хсепіге2 = 600 

40 ycentreX = 500 

50 radius% = 400 

60 sides% = 104 

70 oldsin = 0 

80 oldcos = radius% 

90 FOR 17 = 0 TO sides DIV 8 

100 sin = radiusZ*SINC2xPIX*IZ/sidesA) 

110 cos = radiusZ*COS(2*PI*IZ/sidesZ) 

120 MOVE xcentre%+oldcos,ycentreZ%t+oldsin 
130 DRAW xcentre%+cos,ycentreZ+sin 

140 MOVE xcentreZ%+toldsin,ycentreZ%+oldcos 
150 DRAW xcentreZ%+sin,ycentreZ+cos 

160 MOVE xcentreX-oldsin,ycentreZ-*oldcos 
170 DRAW xcentreX-sin,ycentre+cos 

180 MOVE xcentreX-oldcos,ycentreZ-*oldsin 
190 DRAW xcentreZ-cos,ycentre/-*sin 
200 MOVE xcentreX-oldcos,ycentreZ-oldsin 
210 DRAW xcentreZ-cos,ycentreAX-sin 
220 MOVE xcentreX-oldsin,ycentreZ-oldcos 
230 DRAW хсепігей-сіп,усепегей-со5 
240 MOVE xcentreX*oldsin,ycentreZ-oldcos 
250 DRAW xcentreZ-*sin,ycentreZ-cos 
260 MOVE xcentreZ%+oldcos,ycentreZ-oldsin 
270 DRAW xcentre%+cos,ycentreZ-sin 
280 oldsin = sin 
290 oldcos = cos 
300 NEXT 
310 PRINT 
320 T = TIME/100 
330 PRINT "time taken = ";T;" seconds" 
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9.5 Hybrid Method 


This method combines the previous three methods. The following 
example program will plot a 32-sided circle in 0.7 seconds (under 
0.6 seconds if written in assembly language). Eight plots are 
performed at a time, reducing the size of the main loop to 4. Even 
so, only 1 sine and 1 cosine are calculated, and these are of course 
outside the main loop. 


10 ТІМЕ-0 

20 MODE O 

30 xcentreX = 600 
40 ycentreX = 500 
50 radius% = 400 
60 sides% = 32 

70 oldsin - 0 

80 oldcos = 1 

90 sinc = SIN(2*PI/sides%) 
100 cinc = COS(2xPI/sidesX) 
110 FOR IX = 0 TO sides% DIV 8 


120 newsin = oldsinxcinc + oldcosxsinc 
130 newcos = oldcosxcinc - oldsinxsinc 
140 sinold = radius#xoldsin 
150 sinnew = radius%*newsin 
160 cosold = radius%*oldcos 
170 cosnew = radius%*newcos 
180 X1% = xcentre/-*cosold 
190 X2% = xcentre%ž+sinold 
200 X3% = xcentre%-cosold 
210 X4% = xcentre%-sinold 
220 X5% = xcentre%+cosnew 
230 X6% = xcentreZ+sinnew 
240 X7% = xcentre%-cosnew 
250 X8% = xcentreZ—-sinnew 
260 Y1% = ycentre%+cosold 
270 Y2% = ycentre%ž+sinold 
280 Y3% = ycentre%-cosold 
290 Y4% = ycentre%-sinold 
300 Y5% = ycentre%+cosnew 
310 Y6% = ycentreZ+sinnew 
320 Y7% = ycentre%-cosnew 
330 Y8% = ycentre%-sinnew 


340 MOVE X1%,Y2% 
350 DRAW X5%,Y6% 
360 MOVE X2%,Y1% 
370 DRAW X6%,Y5% 
380 MOVE X4%,Y1% 
390 DRAW X8%,Y5% 
400 MOVE X3%,Y2% 
410 DRAW X7%,Y6% 
420 MOVE X3%,Y4% 
430 DRAW X7%,Y8% 
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440 MOVE X4%,Y3% 
450 DRAW X8%,Y7% 
460 MOVE X2%,Y3% 
470 DRAW X6%,Y7% 
480 MOVE X1%,Y4% 
490 DRAW X5%,Y8% 


500 oldsin = newsin 
510 oldcos = newcos 
520 NEXT 

530 T = TIME/100 

540 PRINT 


550 PRINT "time taken = ";T;" seconds" 
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10 LARGE MACHINE CODE 
PROGRAMS 


А machine code program can occupy all the memory space 
between PAGE and HIMEM, and in some cases other areas below 
PAGE as well. In a typical game program using one of the 20K 
graphics modes, if the game is disk based, the machine code will 
occupy from &1900 to &2FFF. This is a little under 6K. Although 
this size can often be increased by pinching memory areas from 
&400 to &6FF, and from &A00 to &CFF, the overall machine code 
size is still small and code must be economic. 


However, even though this is undoubtedly a problem, it is not the 
only one or even the main one. The reason for this is that the BBC 
BASIC/assembler requires the source code to be present in memory 
for the assembly step. Thus the amount of machine code that can 
be written in a single go is considerably less than the amount of 
memory available. Since the technique described in this book is 
particularly valuable to those wishing to write large machine code 
programs, it is pertinent to consider the various ways in which this 
problem may be overcome. 


All the methods require that the overall program be subdivided 
into a number of smaller modules each of which is small enough to 
assemble. From each of these smaller modules, the assembled 
machine code is extracted, whilst the source code itself is kept 
quite separately in case further amendment is required. There are 
then two choices. Either the individual machine code extracts are 
merged together to produce a single, large machine code program. 
This is the case if the machine code is to be blown on ROM. 
Alternatively, a BASIC program is written which loads the main 
machine code module and passes control to the appropriate entry 
point address. This is the more usual method for tape/disk based 
games, where the initial BASIC program usually incorporates the 
manufacturer's logo. 


The major problem in all this is addressing. Usually a piece of 
machine code assembled in one area of memory has to be moved to 
another area of memory. This is called ‘relocation’. But 6502 
machine code is not usually relocatable in this way and special 
steps have to be taken to overcome this fact. The problem 
subdivides into two distinctly different relocation problems. 


Firstly there is the problem of relocation of an individual module. 
The problems associated with this are called intra-module 
relocation problems. 
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Secondly, and much tougher, is the problem of cross-references 
between individual modules, giving rise to inter-module relocation 
problems. If Module-A calls a subroutine in Module-B, then 
whenever Module-B is moved, the address of that subroutine 
changes and во must all references to it in Module-A. A similar 
problem arises with any data fields which are referenced in more 
than one module. 


10.1 BASIC 2 Relocation 


Those in possession of a BASIC 2 ROM have a tailor-made way of 
solving intra-module relocation problems in the form of an 
extended range of OPT codes. When 4 is added to the usual values 
for ОРТ, the assembler will assemble the program in an area of 
memory defined by the contents of O96, but the code will be 
assembled as though it were located in an area of memory defined 
by the contents of P96. A simple example should illustrate this 
well. The module below contains a simple subroutine to divide a 
number by 2. At line 20, instead of the normal values for pass% of 
0 and 3, 4 has been added to each to specify that relocation is 
wanted. At line 40, O96 has been set appropriately so that the 
machine code will be physically located in the area mc%. At line 
30, P96 has been set to PAGE, so the assembler will assemble the 
machine code as though it were located at PAGE. The BASIC lines 
at 240 to 250 help in saving this machine code. 


10DIM mc% 1000 

20FOR pass% = 4 TO 7 STEP 3 

30P% = PAGE 

400% = mc% 

5ОГОРТ pass% 

60\ routine to divide a number by 2 
70\ on entry the number is in A 

80\ on exit the result is in X and 


90\ remainder 15 іп Ү. 
100.div2 

110 PHA \ save А 

120 LDY #0 \ zeroise remainder 
130 LSRA \ divide А by 2 
140 TAX \ result in X 
150  BCC out \ no remainder 
160 LDY #1 \ set remainder 
170.out 

180 PLA \ restore А 

190  RTS \ bye bye 
200.finish 

2101 

220NEXT pass% 

230PRINT 
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240PRINT " code is at &";"mc%;" and is 6”; 
245PRINT ;~(finish-div2);" bytes Long" 
250PRINT "relocate at &';^div2 

260END 


The result of running this module on a system with disks (PAGE = 
&1900) is shown below. Notice that the machine code physically 
starts at &1B96 and is &B bytes long. The assembly listing shows 
that the machine code has been assembled to run at &1900, 
however. 


1900 ОРТ pass% 

1900 \ routine to divide a number by 2 
1900 \ on entry the number is in A 
1900 \ on exit the result is in X and 
1900 \ remainder 15 іп Ү. 
1900 .div2 

1900 48 PHA \ save A 

1901 А0 00 LDY #0 \ zeroise remainder 
1903 4A LSR A \ divide A by 2 
1904 AA TAX \ result in X 
1905 90 02 BCC out \ no remainder 
1907 AO 01 LDY #1 \ set remainder 
1909 -out 

1909 68 PLA \ restore А 

190A 60 RTS \ bye bye 

190B -finish 


code is at &1B96 and is &B bytes Long 
relocate at &1900 


To save this machine code we must use: 
*SAVE "DIV2" 1B96 +В 1900 


This entire module can be relocated to any valid address that we 
wish, simply by altering line 30 in the original assembly language 
program, re-assembling and saving the machine code once again. 
Notice, though, that only intra-module relocation problems are 
solved in this way. Any other modules that wish to call ‘div2’ will 
have to do so by absolute addressing, such as JSR &1900. If we 
move the start address of 'div2' we create an inter-module 
relocation problem. More will be said of this later, but for the 
moment let us see how BASIC 1 owners can cope with intra- 
module relocation problems. 
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10.2 Intra-Module Relocation Problems 


Many 6502 instructions will relocate quite happily without any 
correction by the programmer. The previous example program, for 
example, has no problems. It is the machine code instructions 
which reference absolute addresses within the domain of the 
overall program which cause problems. The method to be described 
in this section is often used by BASIC 2 owners in preference to the 
method outlined in the previous section. 


The idea behind the method is to identify those bytes within a 
module which are not relocatable. Of course, this can be done 
manually by simple inspection of the program code, but this would 
be prone to error. It is much safer to let the micro identify the bytes 
which will not relocate. АП that is needed is to assemble the 
module at two different addresses and compare the two machine 
code modules byte-for-byte. Those bytes which are different are the 
bytes which will not relocate. А separate section will later suggest 
ways in which the number of un-relocatable bytes can be 
minimised. One of these ways is so important that it needs to be 
mentioned now. It is a good idea to start each module on a page 
boundary. In this way, intra-module relocation problems are 
reduced to high address bytes only. 


Consider the previous example program. By changing it slightly, so 
that the result and remainder are stored in memory areas 
referenced by labels, we introduce some intra-module relocation 
problems. To identify where those problems are, we can assemble 
the program twice at non-overlapping addresses on different page 
boundaries and test which bytes of the machine code have changed. 


10MODE7 

15REM assemble at &2000 

201% = &2000 

30PROCasm 

35PRINT : REM assemble at &2100 

4017 = &2100 

50PROCasm 

55PRINT 

60FOR IX = 0 TO finish-start 

70 IF 2(&2000+1%) = 2(%2100%17) THEN GOTO 90 
80 PRINT "problem at &';~(&2000+IZ%) ; 
85 PRINT " value = &';~?(&2000+1%) 
90 NEXT 

100END 

110DEF PROCasm 

120FOR pass% = 0 TO 3 STEP 3 

130P% = 17 

140L0PT pass% 

150\ routine to divide a number by 2 
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160\ оп entry the number 15 іп А 

170\ оп exit the result is in 'result' and 
180\ remainder 15 іп 'remainder'. 
190.start 

200.result EQUB 0 

210.remainder EQUB 0 


220.div2 

230 РНА \ save A 

240 LDA #0 \ zeroise 

250 STA remainder \ remainder 
260 PLA \ get A again 
270 PHA \ save A again 
280  LSRA \ divide A by 2 
290 STA result \ save result 
300 ВСС out \ no remainder 
310 LDA #1 \ set 

320 STA remainder \ remainder 
330.out 

340 PLA \ restore А 
350.finish 

360  RTS \ bye bye 

3701 

380NEXT pass 

390ENDPROC 


When the program is run, we obtain a list of the problem areas. 
Normally this would be done with OPT set at 2, but in this case 
listings have been obtained for both assemblies so that the reader 
can verify that the results are correct. 


2000 OPT pass% 

2000 \ routine to divide a number by 2 
2000 \ on entry the number is in A 
2000 \ on exit the result is in 'result' and 
2000 \ remainder 15 in 'remainder'. 
2000 .Start 

2000 00 .result EQUB 0 

2001 00 .remainder EQUB 0 

2002 .div2 

2002 48 PHA \ save А 

2005 A9 00 LDA #0 \ zeroise 

2005 8D 01 20 STA remainder А remainder 

2008 68 PLA \ get A again 

2009 48 PHA \ save А again 

200A 4A LSR A \ divide А by 2 
200B 8D OO 20 STA result \ save result 

200E 90 05 BCC out \ no remainder 

2010 A9 01 LDA #1 \ set 

2012 8D 01 20 STA remainder А remainder 

2015 .out 

2015 68 PLA \ restore А 

2016 .finish 

2016 60 RTS \ bye bye 
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2100 ОРТ pass% 


2100 \ routine to divide a number by 2 
2100 \ on entry the number 15 іп A 
2100 N on exit the result is in 'result' and 
2100 \ remainder 15 in 'remainder'. 
2100 .start 

2100 00 .result EQUB 0 

2101 00 .remainder EQUB 0 

2102 .div2 

2102 48 PHA \ save А 

2103 A9 00 LDA 40 N zeroise 

2105 8D 01 21 STA remainder А remainder 

2108 68 PLA \ get A again 

2109 48 PHA \ save A again 

210A 4A LSR A \ divide А by 2 
210B 8D OO 21 STA result \ save result 

210Е 90 05 BCC out \ no remainder 

2110 A9 01 LDA #1 \ set 

2112 8D 01 21 STA remainder \ remainder 

2115 -out 

2115 68 PLA \ restore А 

2116 -finish 

2116 60 RTS \ bye bye 


problem at &2007 value = 820 
problem at &200D value = 820 
problem at 82014 value = 820 


Thus we can relocate this program on any page boundary providing 
we change just 3 bytes. For example, to relocate the machine code 
at &1500, we would simply enter: 


*LOAD "DIV2" 1500 


281507 - 815 
281500 - 815 
281514 = 815 


10.3 Intra-Module General Case 


The previous section illustrated the way in which intra-module 
relocation problems can be identified. Generally, however, there 
will be insufficient memory to assemble the module twice in situ. 
In these cases it is necessary to break the problem into several 
steps. 


Firstly, assemble the module at one page boundary and *SAVE the 
machine code. Secondly, assemble the module at a different page 
boundary and *SAVE the machine code again. Then *LOAD each 
version separately to non-overlapping addresses. Finally compare 
each byte-for-byte and report on differences. 
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step 1. 


The example program is assembled at &2000 and *SAVEd as 
‘DIV2A’. This is done by sandwiching the module between a few 
BASIC instructions. OSCLI (&FFF7) is used for the *SAVE, but it is 
written in such a way that it works on BASIC 1 and 2 ROMs. 


10mod$ = "DIV2A" : begin% = &2000 

30FOR pass% = 0 TO 2 STEP 2 

40P% = begin% 

5OLOPT pass% 

60\ routine to divide a number by 2 

TON on entry the number is in A 

80\ on exit the result 15 in 'result' and 
90\ remainder 15 in 'remainder'. 
100.result EQUB 0 

110.remainder EQUB 0 


120.div2 

130 РНА \ save A 

140 LDA #0 \ zeroise 

150 STA remainder \ remainder 
160 PLA \ get A again 
170 PHA \ save A again 
180 LSR A N divide A by 2 
190 STA result \ save result 
200 ВСС out \ no remainder 
210 LDA #1 \ set 

220 STA remainder \ remainder 
230.out 

240 PLA \ restore А 
250  RTS \ bye bye 

2601 


270NEXT pass% 

280REM «SAVE this machine code 

290DIM XZ 256 

300$X4 = "SAVE " + mod$ + " " + STR$^beginZ + " " + STR$^PZ 
+ " " + STR$^beginZ 

310Y% = XZ DIV 256 

320CALL &FFF7 

330PRINT CHR$133;'"x";$XX : END 


step 2. 


By changing line 10, the module is re-assembled at &2100 and 
*SAVEd as ‘DIV2B’ 


10mod$ = "DIV2B" : begin% = 82100 

30FOR pass% = O TO 2 STEP 2 

40P% = begin% 

5ОГОРТ pass% 

60\ routine to divide a number by 2 

70\ on entry the number is in A 

80\ on exit the result is in 'result' and 
90\ remainder 15 in 'remainder'. 
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100.result EQUB 0 
110.remainder EQUB 0 


120.div2 

130 PHA \ save A 

140 LDA #0 \ zeroise 

150 STA remainder \ remainder 
160 PLA \ get A again 
170 PHA \ save A again 
180 LSR A N divide A by 2 
190 STA result \ save result 
200 ВСС out \ no remainder 
210 LDA #1 \ set 

220 STA remainder \ remainder 
230.out 

240 PLA \ restore А 
250  RTS \ bye bye 

2601 


270NEXT pass% 

280REM «SAVE this machine code 

290DIM XZ 256 

300$X4 = "SAVE " + mod$ + " " + STR$^beginZ + " " + STR$ PZ 
+ " " + STR$^beginZ 

310Y% = XX DIV 256 

320CALL &FFF7 

330PRINT CHR$133;'"x";$XX : END 


step 3. 


The following program, called ‘INTRA, is a general purpose 
program for comparing machine code modules up to &1000 bytes 
in length. It asks the user to provide the first part of the module 
name (in this case ‘DIV2’). It then appends an ‘A’ to this name апа 
loads the machine code module into an area, load96. It then 
appends a “В” to the name and loads this machine code module at 
load%  &1000. Finally it compares the two modules byte-for-byte 
and reports all differences. 


10MODE7 

20REM set aside an area to «LOAD up to 81000 bytes twice 
30DIM Load% &2000 

4OPRINT TAB(12,2) ;CHR$141; CHR$131;" INTRA" 

5OPRINT TAB(12,3) ;CHR$141; CHR$131 ;" INTRA" 

60PRINT TAB(O,6);CHR$134;"enter module name "; 
TOINPUT "> " mod$ 

80A$ = "LOAD " + mod$ + "A " + STR$ Load 

90B$ = "LOAD " + mod$ + "B " + ТКФ” (LoadZ+&1000) 
100PROCoscli (A$) 

110REM get size of module loaded 

120size% = 256x?&2F9 + 2&2F8 

130PROCoscL i (B$) 

140count% = 0 

150VDU2 : REM turn print оп 

160PRINT "Intra-module problem addresses" 
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170PRINT 

180FOR IX = 0 TO sizeX-1 

190 IF ?(LoadZ+IZ) <> ?(LoadZ+IZ+&1000) 
THEN PRINT "&";^IZ;" from start" : 
count% = count% + 1 

200 NEXT 

210PRINT 

220PRINT "total problems = ";count% 

230VDU3 : REM turn printer off 

240END 

250REM ----------------------- 

260REM routine to *LOAD the machine code 

270REM --------------------- 

280DEFPROCoscLi (C$) 

290PRINT 

ЗООРВІМТ CHR$132;"»";C$ 

310DIM XZ 256 

320$X4 = C$ 

330Y% = XX DIV 256 

340CALL &FFF7 

350ENDPROC 


When the program is run on our example modules, it produces the 
following results: 


Intra-module problem addresses 


&7 from start 
&D from start 
814 from start 


total problems = 3 


These are exactly the results that we expected. By inspecting the 
assembler listing for this module, we can determine what value 
with respect to the origin page of the module must be inserted into 
each of these bytes in order to make the module relocatable. In this 
example the values are all the same as the origin page. If this is the 
main module, these fixes must be coded into the initial BASIC 
program. Otherwise they are coded into the main module itself. 
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10.4 Minimising Intra-Module Problems 


Clearly, the less intra-module relocation problems that are present 
in the machine code, the easier it becomes to amend and relocate 
that code. A few simple rules can minimise these problems: 


а) Start each module on a page boundary. 


b) Avoid JMP instructions whenever possible. Use branch 
instructions instead as these are always relocatable. For 
example: 


CLC 
BCC Label 


instead of 


JMP Label 


с) Itis frequently necessary to pass parameters from one 
subroutine to another. Parameters that are passed via zero page 
locations, the resident integer variables, the 6502 registers or 
the 6502 stack cause no problems. However, using data fields 
within the domain of the module will cause problems. 


d) Document and test all code thoroughly before saving the 
machine code. This is an advisable practice, even for single- 
module programs. For multi-module programs it is more-or-less 
obligatory because it is so much more difficult to amend 
program code. It is frequently possible to test individual 
subroutines by driving data through them from BASIC. The 
example programs in this book do this quite a lot. 


10.5 Inter-Module Relocation Problems 


All that has been achieved so far is that individual modules сап be 
made relocatable. The harder problem has now to be addressed; 
how can one module communicate with another module? 


In the discussion that follows, we will assume that only two 
modules are needed for the application. When more modules are 
required, the same principles apply. The principles are, however, 
easier to explain with just two modules. 


Once again, the problem becomes easier to manage when we 
subdivide it into its component parts. The inter-module 
communications problem can be broken into code and data 
problems. In the case of code problems we are concerned with how 
to JMP or JSR from one module to another. This problem becomes 
really simple if we specify certain design constraints on the two 
modules. For example, let us specify that there will be a main 
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module, called ‘MAIN’ and a module of subroutines, called ‘SUBS’. 
In this scheme, MAIN will be located in a memory area that will 
allow it to co-exist with the initial BASIC program. SUBS will 
overlay the initial BASIC program at PAGE and must be loaded by 
MAIN. Let us further specify that: 


MAIN can call routines in MAIN 
MAIN can call routines in SUBS 
SUBS can call routines in SUBS BUT not in MAIN 


Already we have simplified the problem considerably. The key 
steps can now be taken. Let us assume that SUBS contains three 
subroutines (sub1, sub2, sub3) and that these are each called in 
several places from within MAIN. If we use say &70 to specify 
which subroutine we want, then there need only be one entry point 
in SUBS. 


In this example code, the only entry point into SUBS is at the start, 
at 'subs'. The first piece of code resolves which subroutine is 
actually required. In practice, this section of code would probably 
also save the register set. 


10REM SUBS 

15DIM mc% 81000 

20FOR pass% = O TO 2 STEP 2 
30P% = mcZ 

40LOPT pass% 


50.subs 
60 LDA &70 \ get which sub 
70 CMP #1 \ is it sub1 ? 
80 BEQ sub1 \ yes - sub1 
90 CMP #2 \ 15 it sub2 ? 
100 BEQ sub2 \ yes - sub2 
110 CMP #3 \ is it sub3 ? 
120 BEQ sub3 \ yes - sub3 
130 ВЕК \ mistake 
140. sub1 
150 JMP dosub1 \ DO sub1 
160. sub2 
170 JMP dosub2 \ DO sub2 
180.sub3 
190 JMP dosub3 \ DO sub3 
200. dosub1 

: assembler 

: code 


9001:МЕХТ pass% 
910END 


Let us assume that our first guess is that MAIN will be located at 
&2800. We know that MAIN will have to load SUBS at PAGE and 
also handle all its intra-module relocation fixes. We can also design 
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MAIN во that it only calls SUBS іп one place. The entry point into 
this program is at ‘main’. 


10REM MAIN 

20FOR passX = 0 TO 2 STEP 2 

30P% = 82800 

4OLOPT pass% 

50.main 

60 LDA #&83 \ get PAGE 
70 JSR &FFF4 


80  STY callsub+2 \ initialise callsub 
90  STY osblok+2 \ апа osfile parm block 
100 LDX #osblok MOD 256 \ get Lo-byte 
110 ІҮ #osblok DIV 256 N get hi-byte 
120 LDA #&ЕЕ \ LOAD SUBS 
130 JMP (82122 \ via OSFILE 
140.back \ return here 

: code 

: to 

: fix 

: SUBS 

: intra-module 

: relocation 

: problems 

: code 
400 LDA #1 \ call sub1 
410 5ТА 870 


420 JSR callsub 


: code 


500 LDA #2 \ call sub2 
510 5ТА 870 
520 JSR callsub 
: etc. 
: etc. 
700.callsub 
710 JMP 80000 \ overwritten by Line 80 
800.osblok \ OSFILE parm block 
810 ЕОМ osname \ point to file name 
820  EQUD 0 \ Load address (see Line 90) 
850  EQUD back \ xqt address (fudge it) 
840 EQUD 0 
850 EQUD 0 
860.osname 
870 EQUS "SUBS" \ file name 
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880 EQUB &D \ CR 
9001:МЕХТ pass% 
910END 


We have now solved all the inter-module code problems. MAIN 
can freely access all the subroutines in SUBS. 


Data problems are even simpler to handle. We simply make sure 
that there are no data problems. Firstly, all data fields are defined 
in MAIN except for any specific work areas needed by SUBS but of 
which MAIN needs no knowledge. Secondly, only a routine in 
MAIN is allowed to access a data field in MAIN. Thirdly, 
parameters passed between MAIN and SUBS and vice-versa are 
always in zero page memory, resident integer variables, registers or 
the processor stack. An example should clarify matters. 


Let us assume that ‘sub1’ is a routine in SUBS that updates the 
player’s score. The data field ‘score’ is of course in MAIN. Let us 
also assume that ‘score’ is a 16 bit integer, and that &80 and &81 are 
used to communicate the score between MAIN and SUBS. Then the 
code in MAIN might well be: 


LDA #1 \ specify sub1 
STA &70 \ 

LDA score N get lo-byte 
STA &80 \ save it 

LDA score+1 \ get hi-byte 
STA &81 \ save it 

JSR callsub \ get it updated 
LDA &80 \ get Lo-byte 
STA score N put it back 
LDA &81 \ get hi-byte 
STA score*1 \ put it back 


In this way, ‘sub1’ has no knowledge of the field ‘score’ but can 
still function as a subroutine. 


By defining design constraints on the individual modules, we have 
ensured that there are no inter-module relocation problems. The 
intra-module relocation problems are identified by program. 


In practice, each application has to be treated on its own merits. 
Sometimes more modules are needed. In any event, an overall 
strategy such as outlined here will greatly simplify relocation 
problems. The strategy will ensure that the machine code programs 
will work correctly on both tape and disk based systems, providing 
all references to relocation problem fixes are relative to PAGE. 
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10.6 Initial BASIC Program 


Let us assume that we have two modules, MAIN and SUBS 
designed along the lines suggested before. Let us assume that we 
have now discovered that MAIN should be loaded at PAGE+&800 
and SUBS should be loaded at PAGE as was always intended. 


We will assume also that ‘INTRA’ identified the following 
intramodule problems :— 


MAIN (originM = page at which MAIN starts) 


&15 from start should be originM 
878 from start should be originM+2 
&FC from start should be originM+1 


SUBS CoriginS = page at which SUBS starts) 
885 from start should be originS+1 


The fix for SUBS has to be coded into MAIN. The fixes for MAIN 
are coded into the initial BASIC program which would look 
something like this: 


10MODE7 

20A$ = "LOAD MAIN " + STR$~ (PAGE+&800) 
30DIM XZ 256 

40$X% = AS 

50Y% = X% DIV 256 

60CALL &FFF7 

70? (PAGE+&815) = (PAGE DIV 256) + 8 
80? (PAGE+&874) = (PAGE DIV 256) + 10 
90? (PAGE+&8FC) = (PAGE DIV 256) + 9 
100CALL PAGE+&800 

110END 
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THE ADVANCED BASIC ROM USER GUIDE 
FOR THE BBC MICROCOMPUTER 


This book delves deep into the BBC microcomputer BASIC 1 and BASIC 2 
ROMS and comes up with 69 useful subroutines that can be called from an 
assembly language program. The routines cover: 


* 32 bit integer arithmetic 

* floating point arithmetic 

* maths. functions such as sine, cosine, log, square root 

* data conversions 

* random numbers 

The author has programmed commercially for 18 years on a wide variety of 
computers including in more recent years the BBC microcomputer. The book 


attempts to fill some of the important gaps in the microcomputer literature 
and covers in addition: 


* making trigonometry faster 

* writing large, relocatable, machine code programs 

* useful pseudo-directives 

There are many program examples in this book and much more besides. 


The serious programmer of the BBC micro will find this book to be a 
valuable aid. 


Published by 
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153-154 East Road, Cambridge, England. 


